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ABSTRACT 

Aims. Our goal is to study the different morphologies associated to the interaction of stellar winds of AGB stars and red supergiants 
with the interstellar medium (ISM) to follow the fate of the circumstellar matter injected into the interstellar medium. 
Methods. Far-infrared Herschel/PACS images at 70 and 160 /um of a sample of 78 Galactic evolved stars are used to study the (dust) 
emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled 
gas and dust in wind-ISM interactions are used to compare with the observations. 

Results. Four distinct classes of wind-ISM interaction (i.e. "fermata", "eyes", "irregular", and "rings") are identified and basic 
parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. The 
total dust and gas mass inferred from the observed infrared emission is similar to the stellar mass loss over a period of a few thousand 
years, while in most cases it is less than the total ISM mass potentialy swept-up by the wind-ISM interaction. De-projected stand-off 
distances (R ) - defined as the distance between the central star and the nearest point of the interaction region - of the detected bow 
shocks ("fermata" and "eyes") are derived from the PACS images and compared to previous results, model predictions and the sim- 
ulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used 
together to independently derive the local ISM density. 

Conclusions. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar prop- 
erties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection 
of a bow shock - for the distance limited sample - appears to be governed by its physical size as set roughly by the stand-off distance. 
In particular the star's peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or 
detached rings. In most cases the derived ISM densities concur with those typical of the warm neutral and ionised gas in the Galaxy, 
though some cases point towards the presence of cold diffuse clouds. Tentatively, the "eyes" class objects are associated to (visual) 
binaries, while the "rings " generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal 
pulse. 

Key words. Stars: AGB and post- AGB - Stars: supergiants - Stars: mass-loss - Stars: winds, outflows - Stars: circumstellar matter - 
Infrared: ISM - Hydrodynamics 



1. Introduction 

The chemical enrichment of the Universe plays an important 
role in explaining stellar and galaxy evolution. Apart from su- 
pernovae, asymptotic giant branch (AGB) stars and supergiants 
play a crucial role in supplying heavy elements into the inter- 
stellar medium (ISM). It has been well established that evolved 



to: 



N.L.J. 



Cox, 



e-mail: 



Send offprint requests 
nick . cox@ster . kuleuven . be 

* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

** presently at: Department of Physics and Astronomy, University of 
Denver 



stars lose material during their AGB and supergiant phase via 
a dusty wind, and this wind material will eventually be dissi- 
pated into the ISM (e.g. Tielens et al. 2005 ). The presence of 
extended envelopes around AGB stars and supergiants is a di- 
rect result of their stellar mass loss. The properties of these (de- 
tached) shells are directly affected by the mass-loss history of 
low- to intermediate-mass stars and thus contain information on 
this late stage of stellar evolution. 

Both the chemical composition and physical conditions of 
the stellar-wind material returned to the ISM can be altered 
by several processes, such as mixing of nucleosynthesis prod- 
ucts into the envelope, dust formation ( |Ferrarotti""& Gail 2006), 
shock-induced chemistry in the inner wind envelope ( Duari et al. 
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1999]), as well as photo-dissociation due to interstellar cosmic 
rays and UV photons penetrating the outer envelope ( Willa cy &| 
Millar 1997; Dec in et"aL]|2010| ). There is a marked difference 
between the dust and gas forme d in circums tellar environments 
and that present in the ISM (e.g. Jones |2Q0T ). To reconcile these 
differences in the composition of circumstellar and interstellar 
dust additional processes are expected to play an important role. 
One such mechanism involves processing and alteration of dust 
in shocks that can occur when the stellar wind material interacts 
with the local ISM. Shocks are not only important in that they 
can alter the composition of the dusty wind but they also gener- 
ate turbu lence and acoustic noise, which affect the evo lution of 
the ISM ( [Cox & Smith|1974t|McKee & Ostriker|1977fr and are 
believed to be of importance to the initial phases of star forma- 
tion ( |McKee & Cowie|1975[|Spitzer|1982| ). 

IRAS observations already hinted at the existence of low 
temperature circumstellar dust ( |van der Veen & Habing 1988]). 
Clear observational evidence for detached shells around evolved 



dergoes additional processing in wind-ISM shocks and that such 
processing is perhaps more common than previously envisioned. 



stars was found by e.g. |01ofsson et al.| ( |1988| ); |Stencel et al 
1988); [Hawkins| ( p~990| ). Following these discoveries, infrared 



observations with IRAS and ISO revealed a growing number of 
AGB stars with extend ed and sometimes clearly detached ther- 
mal dust emission (e.g. Young et al. 1993b a; Waters et al.|1994| 



|Izumiura et al.|199 6; Hashimoto & Izumiura 1998). At the same 
time, CO radio line emission revealed a number of objects with 
detached gas shells (Olofsson 1996, and references therein). A 
number of detached- shell objects was a l so imaged in the opti- 
cal (e.g. |Gonzalez Delgado et al.||2001| |Maercker etaLl|2010 



Olofsso n et al.| |2010). Interferometric radio maps uncovered 



their detailed structure, which turned out to be remarka ble spher- 
ical thick CO-line-emitting, geometrically-thin shells (Lindqvist 
|et al. || 1999} |01ofsson et al.|[2000| >. These observations showed 
that possibly short phases of intense mass-loss rate are required 
to form such large detached (spher ical) molecular shells of swept 
up wind and interstellar mat erial ([Rowan-Robinson et al.|1 986; 
jZijlstra & Weinberger|2002| jLibert et al.||2007| ). The enhanced 
mas s-loss rate events could be induced by a therma l pulse (van 
|der Veen & Habing|1988||Vassiliadis & Wood|1993l ) or alterna- 
tively, the He core-flash at the end of the RGB phase ( [Dominy 
[1984]). 

IRAS also observed for the first time the infrared signatures 
of bow shocks associated to e.g. ultra-co mpact H n regions ( |van 
Bur en et al.|1 990; Ma c Low et al.|199T ) a nd runa way OB stars 



( |Noriega-Crespo et al.|1997| ). |Stencel et al.| ( |1988| ) presented the 
first IRAS detection of a bow shock ar ound a red su pergiant, 
a Ori These data were re-processed by |Nor iega-Cre spo et al. 
(1997 ) which revealed a detached bow shock at ~6' as well as 
a linear bar at 9'. Only recently, infrared space telescopes such 
as Spitzer and AKARI reveal much more detail on the infrared 
signatures of bow shocks around evolved stars, such as R Hya 



(Uet a et al.||2006b a Ori flUeta et al.||2008| ), R CasjHashimoto 
& Izumiura|1998[ |Ueta et al.|2010| ), and U Hya | |Izumiura etaL 
201 1) . The Herschel Space Observatory now reveals these in- 



frared shells and bow shock regions at the best spatial resolu- 
tion ever. Detached infr ared shells around TT Cyg , UAnt, and 
AQ And are discussed in Kerschbaum et al.| ( |2QlU] ) and on more 
objects in [Kerschbau m et al.| (|2011| ). Bow shocks are reported 
for CWLeo ( |Ladjal et al.||2Q10| ), XHer and TXPsc ( |Jorissen| 
et al.||2011| ), and a Ori (Decin et al. in preparation). For o Cet 



the inner part of the stellar win d bubble bounded a nd formed by 
the termination shock is seen (Mayer et al. 2011). For some of 
these infrared bow shocks, counterparts are detected in UV emis- 
sion (c> Cet: |Martin et al.|2007| CW Leo: |Sahai & Chronopoulos 
|2010| ). These observations indicate that the wind material un- 





Fig. 1. Interaction type "fermata" (Class I). PACS 70 /mi (left) 
and 160 /mi (right). The white horizontal bar in each panel is 
r in length, annotated with the corresponding physical size. All 
panels have north up and east to the left. The black line indicates 
the direction of the space velocity of the star (adopting a scale 
such that 1 kms" 1 corresponds to l"on the image). Note: R Scl 
also has an inner spherical shell (not visible here). One needs to 



multiply flux values by 4.25 
MJy sr" 2 . 



10 4 to convert from Jy arcsec 2 to 
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Fig.l. (continued). Interaction type "fermata" (Class I). PACS Fi § L (continued). Interaction type "fermata" (Class I). PACS 

70 m (left) and 160 m (right). For CW Leo see also Ladjal 70 ^ m ( left ) and 160 ^ m ( n S ht )- Radial P rofiles (azimuthally 

et al. (2010) and for a Ori see Decin et al. (in preparation). The averaged over two opening angle ranges: 85-130° and 270-290°) 

bow shock of *Cyg is at several arcminutes to the south (not are shown for 70 A™ ima g e of W Aql (bottom right). Note the 

readily discernible in the image). presence of turbulent instabilities in y. Cep. 



This paper presents the properties - such as morphology, 
size, and brightness - of detached or spatially extended far- 
infrared emission associated to bow shock interaction regions 
and detached shells around a large sample of AGB stars and red 
supergiants. The paper is structured as follows. First we discuss 
the observations, the data processing and map-reconstruction 



scheme, and present the infrared images in Sect. [2] In Sect. 3.1 
we introduce a simple morphological classification system dif- 
ferentiating observed shapes, in particular we ma ke a dist incti on 
between bow shocks and detached rings. In Sects. |3.2|and|3.3|we 



present radial profiles as well as infrared flux measurements of 
the extended emission. In order to make a qualitative and quan- 
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Fig. 1. (continued). Interaction type "fermata" (Class I). PACS 
70 yum (left) and 160 /mi (right), o Cet, R Leo, RT Vir, and X Her 
reveal RT and/or KH instabilities. For o Cet see also Mayer et al. 
(2011). The sources visible eastward of XHer are background 
galaxies (see also Jorissen et al. 2011). 



Fig. 1. (continued). Interaction type "fermata" (Class I). PACS 
70 jim (left) and 160 /zm (right). V1943Sgr, TMIrc, TXPsc, 
and RCas reveal RT and/or KH instabilities. Note: TXPsc also 
has an inner spherical shell (not visible here). The sources visible 
eastward of TX Psc are background galaxies (see also Jorissen et 
al. 2011). 



titative comparison of the observed bow shocks and detached 
rings we review first the basic physics and relevant stellar and 
interstellar parameters underlying the formation of a bow shock 
or detached shell (Sect.|4|. In Sect. [5] we discuss and collect the 
relevant parameters in order to predict, for example, the size of 
the bow shock. Next, in Sect.[6j we present hydrodynamical sim- 
ulations to elaborate on the origin of the observed morphologies, 
in particular, the features related to (turbulent) instabilities. The 
observations are compared to the theoretical predictions and the 
simulations in Sect. [7] Furthermore, the observations are inter- 
preted in the context of stellar mass-loss and ISM properties as 
well as with respect to the interaction between these two material 
flows. Then we discuss the requirements to observe wind-ISM 
interaction, and thus address the frequent occurrence of bow 
shocks for our entire sample of AGB stars and red supergiants. 
The main results and conclusions are summarised in Sect. 



2. AGB stars and supergiants in the MESS Herschel 
Key Program 



(Mass-loss of Evolved StarS) program 
et al. |2011| ) is a large far-infrared Herschel 
([Pilbratt et al.||2010|T survey aimed at studying the mass-loss 



The MESS 
([Groenewegen 



history of evolved stellar objects, from AGB stars to post- AGB, 
planetary nebulae, luminous blue variables and supernovae. One 
key goal of the MESS program is to resolve the circumstellar 
envelopes around a representative sample of evolved stars (from 
AGB stars to PNe), thereby studying the global evolution of 
mass-loss and circumstellar envelope structure. The details of 
the program, the complete target list, data processing approach, 
and first science results are summarised in |Groenewegen et alT] 
( |20TTT ). 
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Fig. 2. Interaction type "eyes" (Class II). PACS 70 fim (left) and 160 //m (middle). Azimuthally averaged radial profiles are shown 
for 70 and 160 jim (right). Azimuth opening angles adopted for the radial profile are indicated at the top of each radial profile panel. 
Note: U Cam also has an inner spherical shell (not visible here). 



Briefly, infrared scan maps of a large sample of AGB 
stars and supergiants are obtained with the Photodetector Array 
Camera and Spectrometer (PACS; Poglits ch et al. | |20 1 ) using 
the "scan map" observing mode with the "medium" (20" s" 1 ) 
scan speed in the blue (70 /mi) and red (160 /mi) filter setting. To 
create a uniform coverage, avoid striping artefacts, and increase 



redundancy, two observations at orthogonal scan directions are 
systematically concatenated. Scan lengths range from 6 to 34', 
scan leg steps are 77.5' or 155', and repetition factors ranged 
from 2 to 8. 

Data processing was performed applying the standard 
pipeline steps. In particular, a careful correction of the detector 
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Fig. 2. (continued) Interaction type "eyes" (Class II). PACS 70 yum (left) and 160 yum (middle). Azimuthally averaged radial profiles 
are shown for 70 and 160 yum (right). Azimuth opening angles adopted for the radial profile are indicated at the top of each radial 
profile panel. 



signal drifts is critical to reveal faint extended emission struc- 
tures with PACS. After deglitching, we applied two different 
map-making algorithms to assess the quality of the final map, 
and to verify the faint emission associated to the wind-ISM in- 
teraction; these are "PhotProject" with high-pass filter (HIPE 
v.7.0.) and "Scanamorphos" (Roussel 2011). The maps shown 
here are those obtained with the latter method as these are less 
susceptible to the filtering of true extended emission. All frames 
are projected onto an image with a pixel size of 1" and 2" for the 
blue (70 yum) and red (160 yum) bands, respectively. Thus, the fi- 
nal image over- samples the instrumental point- spread functions, 
which have full- width at half-maxima of 5.6" and 1 1.4" at these 
wavelengths, respectively. 

Aperture photometry of the central point-like sources ob- 
served with PACS at 70 and 160 yum shows that both the 
"PhotProject" and "Scanamorphos" maps give flux densities that 
agree within 5%, well within the 15% calibration uncertainties. 
However, the processing with the updated calibration files in 
HIPE 7.0 results in 70 and 160 yum fluxes that are 11% and 
15% lower, respectively, compared to the values presented in 
|Groenewegen et al.| ( [20TT| ) using HIPE 4.4. 

Figures[T]to|4]presents the objects in our sample which reveal 
signs of extended and/or detached dust emission around the cen- 
tral star in either or both the PACS 70 and 160 yum bands. Colour 
versions of these figures are available in the electronic edition 
providing enhanced visibility. The different observed morpholo- 
gies are divided in several classes (see Sect. [3]). Out of 78 AGB 
stars and supergiants imaged with PACS as part of the MESS 
survey, 32 stars show clear evidence of wind-ISM interaction, 
15 stars show detached rings, and 6 show extended irregular 
emission (Tables [T] and [2]). The remaining 30 stars do not show 
evidence for wind-ISM interaction (Table [3]). Tables [T] to [^pro- 
vide basic properties of the observed AGB stars and supergiants; 
IRAS identifier (col. 1), Name (col. 2), Distance (col. 3), Mass- 
loss rate (col. 4), z (col. 5); height above the Galactic plane, hh 
(col. 6); local ISM density, \i (col. 7); proper motion, v* (col. 
8); space velocity, P. A. (col. 9); the proper motion position an- 
gle measured from north to east, i (col. 10); the inclination angle 



of the space motion with respect to the plane of the sky, and v w 
(col. 1 1) the terminal wind velocity. In addition, we also provide 
the predicted stand-off distance, Ro in arcminutes (col. 12) and 
parsec (col. 13), as well as, for Tables 1 and 2, the measured 
stand-off distance Ro in arcminutes (col. 14), parsec (col. 15), 
the position angle 6 (col. 16) and the inferred local ISM density, 
He (col. 17). Columns 18 and 19 give information on the spectral 
type / circumstellar chemistry and binarity, respectively. 

After deconvolution, three stars (R Scl, TX Psc, U Cam) re- 
veal evidence for both a wind-ISM interaction as well as a de- 
tached ring. To measure the ring of TX Psc a de-convolved 2D- 
image was used. Detailed studies of some individual (classes 
of) objects included in the survey here are being presented in 
companion papers. The complex multiple wind-ISM interaction 
shells associated to a Ori are presented, together with comple- 
mentary observations, in Decin et al. (in preparati on). The inter- 
action around X Her and T X Psc is described b y |Jorissen et al. 
pOTT] ) and that of o Cet by |Mayer etaL] pITT] ) . First results on 
the detached ring-like shells around AQ And, TT Cyg and U Ant 
were already discussed by Kerschbau m et al. ( |2010|), while 
CWLeo (IRC+10 216) was discussed previously by Ladjal et al. 
( |2010| ) and |Decin et al.| ( |20TT] ). Also, the objects with both inner 
shells and large detached bow shock regions will be discussed in 
more detail in a forthcoming paper. This notwithstanding, these 
sources are included here to complete the sample of observed 
AGB stars and red supergiants. 
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Fig. 3. Interaction type "rings" (Class III). PACS 70 jim (left) and 160 jim (middle). Azimuthally averaged (360°) radial profiles are 
shown for 70 and/or 160 jim (right). The ring towards S Cas is only faintly visible in the current display, but its presence is verified 
from the radial profiles. For AQ And see also Kerschbaum et al. (2010). 





Distance [arcsec] 




Fig. 3. Interaction type "rings" (Class III). PACS 70 jim (left) and 160 fim (middle). Azimuthally averaged (360°) radial profiles are 
shown for 70 and/or 160 jim (right). For U Ant see also Kerschbaum et al. (2010). The ring towards Y CVn is only faintly visible in 
the current display, but its presence is verified from the previous detection with IRAS ( [Izumiura et al . 1 1 996] ) . 
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Fig. 3. (continued) Interaction type "rings" (Class III). PACS 70 fim (left) and 160 jim (middle). Azimuthally averaged (360°) 
radial profiles are shown for 70 and/or 160 jim (right). For TTCyg see also Kerschbaum et al. (2010). The ring around WHya is 
incomplete and deviates slightly from sphericity, also a jet-like structure is visible in the north-east direction. 
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Fig. 4. Interaction type "irregular" (Class IV). PACS 70 jim (left) and 160 jim (right). These objects all show some evidence of 
irregular, diffuse extended emission. However, their morphology is neither a clear (double) arc nor a ring (possibly simply due to 
observational limitations in sensitivity and/or spatial resolution) and are therefore included in Class IV. 



9 



Cox et al.: Bow shocks and detached shells in the far-infrared 



a 



Ph' 



a 



a 



a 
z 



< 



B a 

r- so 

u u o o u u 



oo OO IT) v ; oo 

t-h O OO *-h 

M M rH On CO 



^ ^ CO LT; LT; 

CO SO CN O CN 

t-H CO CN I t-H <-H 



O O co t-h so 
so in co t-h co xt 

CO CO NO t-h ,-h CN 



, "1 rj rr i P 

i CN CO CN CO 



hh P in 



6 u 6 



u 
u 



oo r- 



a ^ ijh <u 
>n 3 + ^ 

\doo^ + 

>5 *5 SO NO 
<5 ^ 00 00 



p cn oo 
s s s s s 



OO OOoooo OOOOO 



CO 

^ O CN 



CO t-H 



On On 
*-h CN 

NO t 



N N \f NO 
^t -h' t-h l> 

co on m co 



CO NO SO i-H 

co On r-' CN 
no i> cn in 



O CN CO 

in °) <N no 

in M h . 

h m m N 



CN NO On ON CN 
in CN ON t-h' d 
CO CO CO CN xl- 



t> On On CN 
t-h' O ^. <-h 

co cn co in co 



On no in ^t 

oo l> O l> 

,-h' © ,-h' d 

v/-n no in o 

i On no t-h 



O O -hO' 



ro ^£ m 

t-H 1—1 CO 



— On 

CN CN O 



2 o 



: o » no m m 

^CN "■>> ^ ^ ^ G CN 



CN « 

l> HH 

in D 

6 6 



NO ^ 


oo 


NO 


ON 


ON 


CO 






in 


NO 




m 9 


CO 


CN 




CO 


oo 


CN t— h 


t-h CN 


d 








d 






CN 


d 


d 


CO o 


d 


CO 








d CN 




CN 




o 


r- 










i> 


o 




m 










oo 




ON 


"xf 


o 




CN 






oo 


ON 


CO 




CO 


oo 




m 


CN 


co in 




CN 


in 


CN 


,— 1 


OO 






CN 


CN 


CO 




,— 1 


oo 




oo 


CN 


CN NO 
























CO 




















































in 






NO 








CO 




^ 


d 










OO 


o in 




oo 


ON 


CO 


CN 


CO 






in 




CN 


OO On 


NO 


NO 




ON 




in in 


2 o 


P 






P 








P 




O 


O CN 








^ 
















d 








° 


d 


d d 












d d 
























CO 
















i> 


OO 


CN 


m 








in 


m 


o 


NO 




oo 




CO 


m 


CN 




m 


ON 


CO 


CN 


NO 






in 


CN 


NO 


OO \Q 




r- 




oo 


-xt 


no in 


? d 


d 








NO 










d 


d ? 




d 




d 


d 


d t-h' 


CN 




































in 


OO 


ON 


CO 


o 


NO 


CN 


o 


ON 


CN 




o in 




m 


CO 






CO NO 


o ^. 




CN 


CO 






o 




o 


O 


O \f o 


o 




o 


CO 


o 


o o 


VI o 


d 


d 


d 


d 


d 


d 


d 


d 


d 


o 


o o 


d 


d 


d 


d 


d 


d d 


On 




































in ^ 


CN 


r- 


CO 




in 


CO 


NO 






"xT 


o o 


NO 


ON 




oo 


OO 


CO NO 


VI -h' 




i> 


CO 


"xT 




d 


CN 


CN 


d 


d 


t-h' 


d 




d 


CN 




d t-h" 












































































ON 




































CO 












00 






























































o 


q 










U-i 


NO 




© 




q 




q 


"in 


In ^ 


od o. 


o 


^t 




o 


"xt 




ON 


CN 






d *n 






"oo 


in 




d co 


t-H 0O 


i> 






ON 




od 


od 






NO 


CN 0O 


in 






CO 







^ o 

NO -xt 
NO "xf 



u u u o u u u 



t-h q in cn no 
co d d ' — 5 d d d 



oo o 



_ _ _ d d 

cn cn r- oo r- >n on 

t-H ro O ,-H ^H o 

d d d d d d d 



,_; oo |> On 

^ ^ o co C P. in 

ON NO cN SO CO On 

o t-h' d I t-h' t-h' d 



NO t-h NO ON CO 0O 0O 
NO CO t-h CN t-h O O 

d d d d d d d 



CO t-h in 0O CN oo 

in t-h t-h co t-h d d 











o 


^in 


rl "oo 


O 


d 

CN 


d 


^ d 


NO 



^ 00 CN t> On 

O l> ^ On t-h' ^h' on 

>n -xt O i> cn co xl- 

co cn r~ cn co cn co 



i> cn oo q i> <-H 

°°. od od co ■xt so On 
On CN CN CN CN CN CN 



xt t-H ,-H m 



o o o o o o o 



SO CN CN t-h ' 



t-H © 
OO l> ' 

l> ro ( 



< o Z ^ ^ 



2 oo CN £ g oo 

^ ^ ^h 

^ CM CO IJ 5 NO 

+ ro o + ^ 

oo s6 oo r- OO OO 

■xt -xt NO O CN -H 

cn cn t-h in o "xt 

O <— i cn co in in 
OOOOOO 



' CO I +1 <N o © oo 
• in xt rH CN in o 

■ t-H CO l> 0O t-H tpH CN 



- 00 © 

XT) 5h ^ ^ ~ 

t-H t-H NO 0O I 

+1 +1 ^ " +1 11^^ 

NO xt O CO O O ' 
CO CN CO t-h CO t-h ( 







IT) 
CO 


















o 






"oo 


CO 


+1 




+1 


+1 




O 


-xt 


in 


CN 


ON 
CO 




i> 

CN 



on ^t no in 

NO CO t— ' CN 



x t-h in co co -xt 

in so r~~- ■xt 

>o +1 +1 +1 +1 +1 

O On O t-h O SO CO 

CO t-h OO so t*** oo CO 

"xt "xt CO CN CO CO CO 



O £> J 

s ^ Pi 



CO On On 

CN CN CO 

t OO t-H 

O co t-h 

+ + + 

xt t-H 0O 

CN CO xt 

in co xt 

in SO On 

o o o 



O so 
O CN 

CO CO 



X ^ >< 



% 8 g $ 

CO t-h SO so 

O t-h CN 0O 

O O t-H -xt 

■xt NO On On 



XI < w )_!_, -3 

On c3 a <u <q 
> X H =1 W 



cn m m so 

CN O CN qo CN 

l> O OO °5 CN 

cn so cn Y q 

OO in OO On On 

CO t ■ "xt t-h CO 

CN CN CN CN CN 



i u 

: 06 



CO t— H 

O m 

+ + 

oo oo 

co in 

■xt- in 

CO CO 

CN CN 



t3 oo 



5 a «h 

hS « 

>H U >< a 

> Ph > ^ < 



ON xt O , 



CN O Tt g ^ NN ^ 

CN 00 f- S ^ SO S 

no so so 22 !>; » )£ 

+ + + ^ + ^ 

■xt On SO O On CO xt 

t ■ in t-h oo oo co t-h 

co xt xt in co cn co 

co xt O O t** - On On 



10 



73 



O 




Cox et al.: Bow shocks and detached shells in the far-infrared 



11 



Cox et al.: Bow shocks and detached 



shells in the far-infrared 



a 

e 3 



Q B 



B 

z 



u o u u u u u 



O O CN CO CN • 



l l l l l l l 



in in h 

co CN o cN cN 

(N CN O O O 

d> d> d> d> d> d> d> 



r> co co cn 

oo oo CN i-h C""- in On 

o o o o o -h' -h' 



o o o o o o o 



i h \t m ' 



o ^ m . 



o o no £ m 

CN OO CN oo 



p th co oo in 

06 d 06 d ^ in 06 

on o in 00 

CN CN CO ON CN 



O co 00 in On 
^" rn' iO * 06 VO h 

in on on co r- 



r- in co 00 
co on no °i d ON NO 

CO NO 



O O O o o o 



8 vo 

cn 



^ CO 
£ © 2 © 

1 NO 1 ON 



S u 
uouuuuuu 



© no co ^ cn 00 p 
odd h ci in h 



0000 



O O O CN CN ' 
©©©©©©©© 



; NO CN ON CO 

I d cn cn in d 



n S cs NO S Om 
^ ■> ©. CN P CN d 
r-OOOOOO^HOOCN'-H 



no 00 in in h 

CN ^ CN On ^ ^ ^ ^ 

MOoo^inoo^^t 

COCNNOOOCN^H^HCN 



r~; cn ^ in co no 

CN OO NO OO i-H CO NO 

CO^^H^HCOCNCNND 



: tJ- cn cn (-5 
1 in on in co in 



-h o o o 



5 Oi . . M 2 
CN cN 00 CN in 
co f- 1 1 m , -h 



£ CN 



^-cN§ c d^Sg^^^^^^°'Jn < 

NDCNC^-HCN^SdO-HO-H-Hw^-H( 



P CO l> ^ 

O CO NO ^ . _r 

S S S 55 u 

O O O co U 



CO ON ON NO ON 

O O O CN o 

d d d d d 



NO On CN OO 

o -h' ^ ^ d 






00 


NO 






ON 




CO 


00 


NO 




CN 
CN 



CO OO CN NO 



CN CN 



xj- on on ^ in 

^ CO O CO 

^ ^ d d ^ 



> > 
o 

+ - 

in r- 

co ^ 

co 6 



IT) 

00 o 
d d 



CN 
NO CN 

d d 



CO co co ^ ON ^ ^ 



^ 2 H o5 a 

O m m " in no 

VhC2^ od d 



IT) 

CN 
QC 


O 

NO NO 


5 


O CM 




CN 


CN 


NO ^ in 00 
^ CN 1—1 j> 


In 
in 


in -e " js " 
CO ON © 
H W H 


■* ■* -« " -s o-« t"- »s 

incNt^iniicNOO 
co^Hinoovnin^Hco 


* *" 2 

t 00 H © ^ 

CN CN in CN 


781 + 
943 


<n ^> +1+1 +1 

r- 00 00 
r- co r-- no © 
co co fN <N 


+1 +1 +1 +1 „ +1 +1 +1 
H^o^o^Hhin 
cNONoooriONor- 

CO^HCOCOTtCN^j-CN 


+1 +1 +1 <n " +12 

xt-©co©©NOr>© 
r-OONCOCNNONOirj 



u co U u < 53 
< w pi p ^ ^ 



> ts u u 

^^Xc^hScoH 



, 00 r- 



< 



+ + °? + + °? 

00 ON NO ^ 00 ON O 

xl- in ^ r- cn cn in 

CN i-H CN CO O CO CO 

O co in O O 

O O O O O 



in 



r- co 00 

O in in 

00 On i> 

CN NO O 



CO CN 
CN 

00 CO 

r- o 



i n VU U _|_ 1 

t-~ CN "sh NO d c<; c.x; 

cNNOONr-ON^i-inco 
cNcoinooONO^Hco 

^H^H^H^H^HfMCSICN 



CN u u 

o s ffi w u 

O H N ^ 

< o CO p< > 



f- m 

m rt- 00 2 » 

^t- -h Tt a ^ 

+ + + T + 

NO F"- CN O NO 

in on on cn on 



CN CO 

l-H CO 

no in 



CN CO 
CN CN 



12 



Cox et al.: Bow shocks and detached shells in the far-infrared 



X 



as 



0) 0) 0) 0) 



IT) OO OO 

s s s 



13 

co ^o 



600000060660000600600000006 



VI VI o O O VI o O VI o o vi o 



in h o \o h > ^ 
° h «! h d vi 
O VI VI O VI VI CM 



^ooinONOo^ONOi^cooocN^coincN^o ^ 00 ^o co 00 On in 

H O O O ^' O d H H ^ (S O H \D 9 H ^1 H IO CO CO rH 

vi vi w, ^ C;, VI O co O VI VI VI VI 



0<N tN^ " H ^ O V1h OOO 
^HONCNOOCN^H^H^H^HOO^H^HCO 



OO IT) 

© °°. »n l> ^0 

(N OO O 

<-h CN CN 



CO ^O OO On On 

°^ <sj so no c4 -h' ^1 

no »n otN^H-xf 10 00 

OO ^^H^H^H CN 



O co 00 O (\| O tJ- O OO 

Al Al Al CN in co co Al -h ^ -h Al CN CN ^ in Al -h 00 Al CN Al Al in Al Al Al 



•st- m ^o ^o r- 

ir : cn 00* °i i> co 06 
\o ^ ^ r- — cn <n 



00 co o • 

~ — , O O V 

^'ddd^oo^oooooo- 



Ooo^H(N^(N^HOmin 
d ^ d d ^h' d d d d ^ 



in n in S n R^H^oooHir) vo in cn o co cn r-^ 
im^Hi(Ni ■ c — ■ 1 HcfimNmHininHoiom 



in 



co 



On 



o m o 

- CO On 



O 00 O ^ c 
On co O H 
h >o in 0\ 



o <n O on co , , . . 

3 10 in co 00 r-- on °°. ~. ^ 00 — in t-~ in 

*0 (N M rH M CN CO [-» VO T^- 'xj" 



O On O 

^ r- 

l> CN CN 



g £ ffl 

^ < o 



^ co -f lr ' ' 3 ■« o ^ © r- 

OOCO C~- " 1-H ~ ^ rH ff) "fi 00 

2 +1 +1 co M +1 o +1 S S +1 +1 +1 2 £ +1 2 

O CN On +IOcO-xtfOOOO-HOOOCNO 



o 



+1 d c 

On r- ( 
1— 1 CO I 



d °^ 
+ £ 



to 




QC 
ON 


© 

fO 

NO 


ON 

+1 


2 * 


t " 


CO 

5 


O 



00 
l> 



. o o w> 'S O 

a y ~ ? g 



K>><HOS^JJZ 



< 



On OO 

h in h 

CN NO CN 

NO ^O 



ONO^ogcNoS^ONOO^ft. 

HminH^in^ iinH^^^ 
^or-mcN^^m^^o^Co^' 
M >n ^ vd v , 1 - M x 4. <n ^- + + 1 



^CO^^JO^ 

+^ + +??^? + + 

ocNO^^cooN^T^^r^^oN^r^oo^oor--r---^oo 
t^^(N^^r^or^cooNco^^ONi^in^o^r---ON^ocN 
OHOmininNOH^cfiinNHTtoNmmOHHm 
cNco^^^^r^ONOO^HCN-vtininooooooooo^-icoco 

>00000000-H-H-H-H-H-H-H-H-H-H-HCNfN(N(N 



660 



OO o ^o 
in in h 
d d ^h' 



+l +l ._ 

co co On 

in ^j- 

CN 



CO 

>H 00 < 



CN rf 

^ co in 
O co 



3u 
Is 

e3 a 



Cox et al.: Bow shocks and detached shells in the far-infrared 



3. Observational diagnostics of wind-ISM 
interaction regions 

3.1. Morphological classification of the detached far-infrared 
emission 

From the observations (Figs. [T] to [4]) one can immediately recog- 
nise different overall shapes of the detected (detached) extended 
emission. The two most obvious cases are arcs or "fermata" 
and "rings". Two additional distinct morphologies are the dou- 
ble opposing arcs or "eye" and the "irregular" emission. The 
morphological classification is summarised in Tableland indi- 
cated for all sources in Tables Q] and |2 

"Fermata" (Class I) interaction is characterised by a rela- 
tively smooth elliptical arc spanning an azimuthal opening angle 
>120°. This class resembles closest the wind-ISM bow shock 
shape predicted by theoretical models (Sect. [4]). We note that 
there are different cases within this class reminiscent of the dif- 
ferent shapes that occur due to variations in the stellar and inter- 
stellar parameters. This will be discussed in detail in Sect. [6] For 
example, X Pav, EP Aqr and X Her show a peculiar bullet-like 
shape with back flow emission in the wake of the bow shock. 
Others show clear signatures of Kelvin-Helmholtz (KH) and/or 
Rayleigh-Taylor (RT) instabilities (see Sect.|6]). Class 11 ("eyes") 
includes objects with two elliptical non-concentric arcs observed 
at opposing sides of the central source, both have a covering an- 
gle of < 180°. In two cases, VYUMa and UCam, the two arcs 
are connected and there is even tentative evidence for a jet struc- 
ture in the mid plane. Class III consists of spherical structures, 
i.e. circular "rings". This category includes typical detached 
spherical- shell objects, but could also include true rings (i.e. not 
spheres). From the seven, well-studied CO-detected detached 
shells around carbon star s (e.g. UAnt, UCam, TTCyg, RScl, 
S Set, V644Sco, DRSer; |Maercker et al.pOTO] ), only the latter 
two are not included in the MESS survey. All known and ob- 
served detached CO shell objects show a spatially resolved dust 
ring, co-spatial with the gas emission. Larger, and thus presum- 
ably older, rings are also detected around Y CVn and AQ And. 
For these objects no corresponding CO shell has been detected, 
possibly due to the photo-dissociation of CO by the interstel- 
lar radiation field (Libert et al. 2007; Kerschba um et al.||2010| ). 
Sources with diffuse irregular extended emission are classified 
"irregular" (Class IV). All other targets, for which no evidence 
of diffuse extended emission has been observed, are assigned to 
Class X, "non- detection " . 

Interestingly, for several objects, such as RScl, UCam, 
TX Psc we detect small detached rings in addition to the classical 
bow shock region further away from the central source. These 
objects are subsequently assigned to both Class I and Class III 
(Tables [T] and [2]). For these objects it may be the case that we 
observe a young spherical shell (originating from a wind-wind 
interaction due to a recent thermal pulse) expanding in a rela- 
tively low density medium. In this scenario the local environ- 
ment of the star has been blown out by the earlier wind/mass- 
loss which swept out ISM material and created a wind-ISM 
bow shock. At a later stage a thermal pulse produced a den- 
sity / temperature enhancement traced by the infrared emission. 
Alternatively, the small inner ring represents a structure delineat- 
ing the interface between the (back flowing) termination shock 
and the free expansion zone. For other objects, such as CWLeo 
and a Ori, the Herschel/PACS observations reveal irregular mul- 
tiple incomplete shells in the inner regions of the stellar wind 
envelope (e.g. Decin et al.|2011| ). 



Table 4. Morphological classification. 



Class Description 



Shape 



I Fermata a 

II Eye - double non-concentric arcs 

III (circular) Ring 

IV Irregular (diffuse) 



Notes. {a) A "fermata" is a musical sign. 




Fig. 5. Illustration of the extended flux measurements of the ob- 
served bow shock of 6 Aps. Both the elliptical bow shock aper- 
ture and the sky aperture are indicated. 



3.2. Radial distance of arcs and detached shells 

For some objects the extended emission is very faint with a low 
contrast with respect to the sky background. In order to improve 
the measurement of the angular distance/size of arcs and shells 
we have constructed radial profiles for Class II and Class III ob- 
jects. Table [5] gives the radii of the "eyes" (II) and "rings" (III) 
as identified in the azimuthally averaged radial profiles shown 
next to the image panels in Figs. [2] and [3] (adopted azimuth an- 
gles are given at the top of each radial profile panel; for the 
"rings" objects it is alway 360°. For the "eyes" the radii and 
radial profiles of both arcs are given. For a few cases with very 
faint extended emission, the assignment is ambiguous. VPav 
possibly belongs to Class III instead of Class II as the dust emis- 
sion faintly traces a full circle. RTCap is tentatively assigned 
to Class III but its ring is not complete and rather consists of 
two distinct arcs. However, the arcs are concentric, which is not 
the case for true "eye" objects. For SCep the ring is also not 
complete but resembles more a full circle with gaps. The de- 
tached shell around Y CVn is very faint and barely visible at a 
radius of -190", but in this case the detection is supported by 
the radial profile and its e arlier detection with ISO/ISOPHOT 
by |Izumiura et al. (1996]). Tentatively, a faint arc bow shock 
structure ("fermata") extending over -100° (from position an- 
gle -30° to 130°) is present at 9' east of the central star. XTrA 
is also listed in Class III although also here the ring is faint, with 
only a brighter arc to the east. Further observations are necessary 
to resolve these ambiguities. 
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Table 5. Aperture flux of observed bow shocks (Class I and II) and detached rings (Class III). Irregular Class IV is excluded. An 
elliptical annulus was applied in order to measure flux in the non-spherical arcs. The annuli of the detached rings were nearly 
spherical. The shape of the elliptical annulus was identified from the image with the most significant detection and then directly 
applied to the other band. For the fermata the flux is only integrated over a limited azimuthal angle covering the observed emission. 



IRAS id 


Object 


Class 


Radius a 


Dust annulus 


Flux (Jy) 


-^dust+gas 


(10 M ) 












(") 


(pc) 


(") 


70 fim 


160 fim 


70 //m 


loll /im 


(10 Mo) 


/1 n-4 A/r \ 
(10 Mo) 


00248+3518 


AQ And 


I+III 


52 


0.21 


40-62 (circle) 


2.20+0.01 


1.28+0.01 


0.9 


1.1 


1.1 


4.4 


ai i cn . noon 

01159+ /zzo 


S Cas 


III 


50 e 


0.23 


32-64 


2.26+0.01 


1.11+0.02 


1.1 


1.1 


12 


54 


01246-3248 


RScl 


I 

III 


54 j 
14/ 


0.10 
0.03 


51-68 (arc) 


1.27+0.01 


0.49+0.01 


0.2 


0.2 


0.1 


0.3 


02168-0312 


o Cet 


I 


82' 


0.04 


70-150 (arc) 


47.64+0.03 


9.75+0.05 


2.2 


1.0 


0.5 


0.2 


03374+6229 


U Cam 


III 


If 


0.02 




















II 


57/62' 


0.12/0.13 


80-140 


2.87+0.02 


2.44+0.05 


0.6 


1.1 


40 


1110 


03507+1115 


NML Tau 


I 


85 


0.10 


95-130 (arc) 


5.23+0.02 


3.04+0.03 


0.6 


0.8 


3.5 


62 


04459+6804 


ST Cam 


II 


67/84' 


0.14/0.17 


84-122 


1.28+0.02 


0.61+0.03 


0.3 


0.3 


3.1 


2.2 


05028+0106 


W Ori 


I/III 


92 


0.17 


70-120 
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0.3 


8.2 


32 


05418-4628 


W Pic 


I 


34^' 


0.08 


62-90 


0.92+0.01 
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0.1 
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13 


05524+0723 


a Ori 


I 


397;' 


0.25 


510-660 (arc) 


56.68+0.19 


22.64+0.38 


3.8 


3.2 
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06331+3829 


UU Aur 
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82^' 


0.14 


100-140 (arc) 


5.44+0.02 
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0.9 


0.9 


14 


88 


09448+1139 


R Leo 


I 


93^' 
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94-134 
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0.1 
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2.1 


09452+1330 


CW Le</ 
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I 
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0.29 
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0.8 


76 


31 
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42 
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0.4 
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0.12 
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RT Vir 
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I 


96^ 


0.05 


200-245 (arc) 


4.65+0.02 
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0.3 
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12 
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II 
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0.6 


21 


12 
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S Set 
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70 
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(\ c\r\ / 1 \ 

30-90 (circle) 


14.09+0.02 


8.85+0.03 


2.7 


3.7 


13 
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45/75'' 
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0.5 
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I 
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0.6 


0.5 


3.1 


19 
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III 


92 


0.13 
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15 
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11 
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1W 


0.15 
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11.5+0.12 
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4.65+0.01 
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0.7 
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100-160 (arc) 


8.35+0.01 


3.52+0.04 


0.5 


0.5 


2.5 


2.9 



Notes. (a) Radii of rings and fermata derived from the azimuthally averaged radial profiles (Fig. 1) (b) Derived from Eq.hjusing the total integrated 
flux at 70 and 160 //m, respectively, adopting a gas-to-dust ratio of 200 and a dust temperature of 30 K. (c) M IS m = \kt a>ism> with n H taken from 
Eq.[5] {d) Mism, with n H taken from the local densities inferred from the measured stand-off distance and Eq.[2j Tables [l]and|2] (e) Central source is 
offset by 0" Ra, 5" Dec. if) Ring radius from the deconvolved image. ig) Inner and outer (at 160 11m) rings, respectively. Central source is offset 
by 12" Ra, 5" Dec. (z) Radial distances are quoted for both north and south arcs (east- west for W Aql). (7) projected distance A (Fig. [6]). CWLeo 
has been observed two times at 160 fim. 

to the low densities, T^ ust «: r gas . Low velocity shocks on the 
other hand give primarily rise to dust emission in the infrared. 
Despite observational efforts with Spitzer and Herschel, no in- 
frared line emission has been detected in AGB star bow shocks 
yet (Uet a|2011[ Decin et al. 201 1), and we will assume here that 
the observed emission is entirely due to thermal dust emission 
(see also Sect. [4}. 

In Table [5] we give the aperture flux measured in both PACS 
bands for the elliptical bow shocks and circular rings. To deter- 
mine the appropriate aperture we fit the extended dust emission 
with an ellipse. The inner and outer annulus of the aperture are 
subsequently derived from the azimuthally averaged radial pro- 
file. For detached rings the fitted ellipse is circular and the aper- 
ture is taken over the entire annulus. Bow shocks have a nearly 



3.3. Far-infrared dust emission 

Cold dust grains emit strongly at mid- to far-infrared wave- 
lengths. Dust emission can, in a simple approximation, be de- 
scribed as a blackbody modified by the grain emissivity; F v oc 
By - Qy, with Q v « v^, where ft depends on the type of dust con- 
sidered. For example, J3 = 1.1 for amorphous carbon in the ISM 
( Rouleau & Martin| 1 99 1 \ as well as for carbon grains in carbon 



stars (Jura 1986), while astronomical silicates have J3 = 2 ( |Volk 
|& Kwok|1988T 

Only for fast shocks (v S h 0C k > 50 km s" 1 ) high shock temper- 
atures are reached that give rise to UV radiation (e.g. observed 
for the bow shock of CW Leo) and possibly a strong [0 1] 63 yum 
line. But, as the grains are only weakly coupled to the gas, due 
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elliptical shape, and their aperture flux is determined over a lim- 
ited azimuthal angular range of the annulus. An illustrative ex- 
ample of a bow shock aperture measurement is shown for 6 Aps 
in Fig. [5] Table [5] shows that the 70 over 160 /mi flux density 
(F v ) ratio varies between 1.2 and 3.5, corresponding to a dust 
temperature of -30 to 40 K (for J3 = 2). 

The observed infrared flux - for optically thin emission - 
gives a direct measure of the emitting dust mass, and thus - given 
a gas-to-dust ratio - the total mass. Assuming specific dust prop- 
erties one can write ( Li 2005): 



Hn regions due to stellar winds ([van Buren et al.||1990| Mac 



d F i 



d 2 F v 



dust 



k(A)B a (T) k(v)B v (T) 



(1) 



with F v the observed flux (in erg s" 1 cm -2 Hz -1 ), A the wave- 
length (in cm), d the distance (in cm), T& ust the dust temperature 
(in K), By the Planck black body curve in frequency units, and 
andA: y the dust opacity at the observed wavelengths. We adopt a 
dust o pacity of 60 and 10 cm 2 g" 1 at 70 and 160 yum, respectively 



dust opacity ol oil and 11) cm" g at /U and lol) fim, respet 
(from |Li & Draine|200l| who give k(A) = 2.92l0 5 (A/jim) 



P - 2). For the dust temperature we adopt 35 K (as derived 
above), which means the dust is slig htly heated with respect to 
the cirrus dust (7d us t ~15 to 20 K; Boulanger 2001). The es- 
timated total dust and gas masses (adopting an average gas-to- 
dust ratio of 200) emitting in the bow shock region or detached 
shell are summarised in TableB] The total observed masses range 
from ~ 1 x 10" 5 to ~ 6 x M . We note that the dust mass 
inferred from the measured infrared emission is sensitivy to the 
dust temperature, which is not well constraint by the avaialbe 
data, the dust emissivity law (k(A)), which is uncertain by order 
of magnitude and strongly dependent on the chemical composi- 
tion of the dust, as well as the gas-to-dust ratio. 

The total mass of the ambient medium that could be swept 
up by the stellar wind is roughly the volume set by the stand- 
off distance of the bow shock or the radius of the detached ring, 
r, times the mass volume density: M S h e ii,isM = 4/37T r 3 Pism- 
In Table [5] we give M IS m using both the ISM densities, n H , de- 
rived using eq.[5]as well as using n H inferred from the measured 
stand-off distance and Eq. [2] For may cases the observed gas & 
dust mass is less than the potential swept-up ISM mass, but in a 
few cases the derived masses are similar. The potentially swept- 
up ISM mass, Mism ranges from 10" 5 to 10" 1 M Q . In compari- 
son, assuming constant mass-loss rates, the stellar mass loss after 
1000 years amounts too 1 - 50 x 10" 4 M , again similar to the 
observed total mass. This consistency is not surprising, since to 
first order, the total mass accumulated in the bow shock depends 
on its age (i.e. mass of stellar material piled-up) and the ISM 
density (mass of swept-up ambient ISM material). However, for 
bow shocks (spanning a limited azimuthal angle), only a fraction 
- typically less than half - of the stellar wind and ISM material 
is entrained into the bow shock region. Furthermore, both ISM 
and wind material flow from the bow shock apex along the con- 
tact discontinuity to be finally shed in the tails. For detached 
shells, both the ISM and stellar wind mass are presumably pre- 
served, unless either formation or destruction of dust grains and 
molecules alters the gas-to-dust ratio and its chemical composi- 
tion (e.g. photo-destruction of CO by the ISRF or processing of 
grains in the shock). 

4. Interaction between stellar winds and the ISM 

Bow shocks are common in astrophysical contexts and can oc- 
cur where two material flows - of different density, velocity, or 
viscosity - collide. Examples are bow shocks around compact 



Low et al.|1991 Raga et al.|1995|), but also due to wind s in bi- 
nary systems ([Stevens et al.|1992[|Parkin & Pittard|2008|[Pittara^ 
2009 ; van Ma rie et al.|2011a| Cor movement of an object with a 
magnetic field through a medium ( [Baranov et al.||1971 |Cordes 
et al.||1993|), between slow a nd fast winds ( [Steffen et al.||1998| 
Steffen & Schonberner 2000), and between stellar winds and the 
surrounding ISM (Matsud a et al.| |T989 ; Brighenti & D'Ercole 
T995[|Comeron & Kaper|1998| ). 

In this section we will discuss the formation of bow shocks 
due to interaction of a stellar wind with a low density medium. 
First the general case of a star with a stellar wind moving through 
the ISM is considered. We discuss also the special case of a sta- 
tionary star with a stellar wind expanding into either an older 
stellar wind or the interstellar medium. 



4.1. The size and shape of a bow shock due to "wind-ISM" 
interaction 

If a supersonic stellar wind and the ambient medium collide, a 
bow shock interface can be created at a distance from the mov- 
ing star where the ram pressures and momentum fluxes of the 
wind and the ISM balance each other. This is called the contact 
discontinuity. It is a result of interaction between two fluids, but 
shocks are not necessary. In the absence of shocks, there will be 
only collisional heating and no Ha emission. For shocks to oc- 
cur the Mach number M of the fluid velocity v (M = v/v soim d) 
should be higher than unity. I.e. the relative velocity between 
the stellar wind and ISM needs to be higher than the speed of 
sound in the ISM. Note however that the speed of sound in 
the low density isothermal warm neutral medium (WNM) is 
very low (v soim d ~ 1 kms -1 ) and thus even a slow AGB wind 
(v w ~ 10 kms -1 ) will move supersonically. Because the sound 
speed (adiabatic case; y = 5/3) scales with (P/p) l/2 (or T 1/2) ) 
the sound speed will be lower in the diffuse cold neutral medium 
(CNM). In the non-adiabatic regime the ram pressure balance is 
given by p w v 2 = Pismv 2 sm . Assuming that the layers mix and that 
the post-shock cooling is efficient (i.e. instantaneous cooling), 
the thickness of the dense shell is negligible with respect to the 
distance from the star. This is, for example, valid for a slow stel- 
lar wind interacting with a hot low-density medium ( [Borkowski 
et al.|1992| ). In this approximation the stand-off distance, Rq, de- 
fined as the distance between the star and the apex of the contact 
discontinuity or bow shock, is given by (e.g. Baranov et al.|197l{ 



Dyson| 1975[ |Raga et al.| 1995) |Wilkin| 1995) 



' 47T PlSM V 2 



(2) 



with M the rate of mass loss and v w the velocity of the isotropic 
stellar wind (with respect to the rest-frame of the star), pism 
the mass density of the ambient ISM, v* the relative space ve- 
locity of the star with respect to the ISM. The standard bow 
shock morphology consists of a forward shock separating the 
unshocked and shocked ISM, a wind termination shock that sep- 
arates the free- streaming wind from the shocked wind, and be- 
tween them a contact discontin uity separating the s hocked wind 
from the shocked ISM (see e.g. We aver et al . 1 1 977 1 and |Lamers 
& Cassinelli|1999| ). If the cooling of the shocked stellar wind is 
inefficient, a thick hot, low-density gas layer will exist between 
the free-flowing wind and the bow shock. The termination (re- 
verse) shock and the bow shock, delineating the contact disconti- 
nuity, both travel away from this contact discontinuity; the bow 
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Fig. 6. Observable projected outlines of wilkinoids for an Ro of 1 
and star-ISM motion inclinations of i = 0, 25, 50, and 75° (with 
respect to the line of sight). The stellar position is indicated by 
a star. A is the observed (projected) stand-off distance between 
the star and the bow shock apex. B is the observed (projected) 
distance between the star and the bow shock for 6 = 90 °. 6 
for any point x on the projected outline is the angle between the 
star-apex line and the line connecting the star with the point x. 



shock in forward direction and the termination shock in back- 
ward direction with respect to the relative motion of the star in a 
stationary ISM. 

As mentioned above, for efficient cooling, the physical thick- 
ness of the shocked region remains small and thus the two shock 
fronts delineating the contact discontinuity are not resolved. In 
Sect. [6] we present hydrodynamical simulations of wind-ISM in- 
teraction which explore the effect of varying the parameters used 
in Eq.[2](M, v*, nu) as well as other parameters such as the dust- 
to-gas ratio and the temperature of the ISM. These simulations 
also confirm that, despite the formation of turbulent instabilities 
and inefficient cooling in some cases, Eq. [2] gives in general a 
rather accurate prediction for the stand-off distance. 

Adopting appropriate values for the mass-loss rate, stellar 
wind velocity, the star's velocity, and the ISM density (Sect. [5]), 
we can directly predict the stand-off distance, Ro from Eq. 2] 
This predicted value for Rq can then be compared directly to the 
measured Ro obtained by de-projecting the observed minimum 
distance between the star and the bow shock outline. Vice versa, 
we can use the measured Ro and Eq. [2] to derive tzism- This is 
important as the local ISM density is difficult to determine ob- 
servationally (see Sect. [5]). 

For example, let us apply this to a Ori. The star is at a dis- 
tance, d- 197 pc (z = -14 pc), movi ng at an intermediate veloc- 
ity of 28.3 k m s" 1 through the ISM ( |Ueta et al.|2008[|van Marie 
etal.|2011b|). Its mass-loss rate is 3 x 10"° M Q yr" A and the stellar 

(Tabled. From the observed de 



wind velocity, v w = 14.5 km s" 

projected stand-off distance, Ro of 5.0' or 0.3 pc, the local ISM 
density is then found to be 4.2 cm -3 . This estimate is a factor of 
two higher than the density derived from the models of the ISM 
(Eq.[5p2// = 1.7 cm -3 ). The predicted stand-off distance is close 
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Fig. 7. Observable bow shock properties A and B as a function 
of inclination i for a fixed Ro - 1 . 



The 2D- shape of the contact discontinuity (and thus also the 
termination shock and bow shock) can be solved analytically in 
the optically thin approximation ( [Mac Low et al.|[T991[ [Raga 
|etal.|1995||Wilkin|1996| ): 



R(0) = Ro 



V3(l -6 cot 6) 
sin# 



(3) 



with 6 the latitudinal angle from the apex of the bow shock as 
seen from the position of the central star, and Ro, the stand- 
off distance defined in Eq. [2] Alternatively, it can be written as 

y(z)= ^Ro^l-z/Ro (e.g. 



Mac Low et al. 



1991 



Fig.g.Eq.g 



is valid only for the case of a relative star-ISM motion in the 
plane of the sky (i.e. a "side view" of the bow shock or a star- 
ISM motion inclination of i = 0° with respect to the plane of the 
sky). In this case, Ro is directly the angular separation, A, of the 
star to the observed bow shock outline (Fig. [6]). Similarly, when 
there is no inclination^' = 0°), the angular distance between the 
star and the surface of the bow shock parabola at an angle of 
= 90° (perpendicular to the apex direction from th e star) can be 
defined as £, where B = R(90°) = y(0) = R V3 (|Wilkin||l996 



to that obtained fro m a detailed simulation of a Ori (see e.g. |van 
|Marle et al.|2011b| and simulation A in Sect. [6}. 



and Fig. [6]). In this geometry, Ro can also be derived from mea- 
suring the angular distance between the star and the bow shock 
interface at an angle of = 90° with respect to the star-ISM mo- 
tion di rection (B\ as illustrated in Fig. [6]), since R(90°) = y(0) = 

Ro V3 dWilkinj 19961 ). 

In general the relative star-ISM motion will be inclined with 
respect to the plane of the sky, thus i ± 0. In order to quantify the 
effects of an inclined bow shock with respect to the plane of the 
sky, we simulated the appearance of the three dimensional bow 
shock surface. This hollow paraboloid wilkinoid surface is de- 
fined by rotating the two dimensional shape, i.e. the R(#) curve, 
around its axis of symmetry (which is defined by the line con- 
necting the star and the bow shock apex). We numerically con- 
struct a 3D wilkinoid surface. These points are then projected 
onto the plane of the sky under a given inclination angle to com- 
pute the observable outline of the bow shock in the optically 
thin approximation. This outline thus traces the loci of tangen- 
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tial line of sight to the bow shock paraboloid at a given incli- 
nation angle. Four such constructed wilkinoids with an Ro of 1 
and star-ISM motion inclinations with respect to the plane of the 
sky of i = 0, 25, 50, and 70° are shown in Fig.[6]to illustrate the 
change in shape with increasing inclination angle. These simu- 
lations show indeed that for an increasingly higher inclination 
i (corresponding to a smaller viewing angle with respect to the 
line-of-sight), the projected size of the bow shock outline be- 
comes increasingly larger, which is qualitatively consistent with 



earlier bow shock simulations for compact Hn regions (Mac 
Low et al.|1991| ). As noted above, in Fig. [6] we also indicate the 



two observable quantities, A and B. A is the projected - minimal 
- distance between the star and the bow shock outline in the di- 
rection of relative motion, which can be measured directly from 
the observed image. B is the distance perpendicular to A, i.e. at 
= 90°, from the star to the bow shock outline. It is important 
to note here that due to the inclination of the bow shock surface 
with respect to the line-of-sight, a "pseudo" apex of the observ- 
able outline appears where the line-of-sight becomes tangent to 
the rotated 3D-Wilkin paraboloid. The observable A is in fact 
the distance from the star to this "pseudo"-apex projected onto 
the plane of the sky (i.e. not to the real apex of the wilkinoid). 
At high inclination, but less than 90°, both the "pseudo"-apex 
and true apex of the observable outline would appear - in the 
model - as two brightness enhancements. We refer the reader to 
Figs. 5 and 11 in Mac Low et al.|([T991) for an excellent illus- 



tration of this effect. However, due to superposition of the true 
apex and the very bright central star the latter will not be de- 
tectable for our stars. In addition, the faint "pseudo"-apex will 
also be difficult to observe. For i = 90° there is no tangent to 
the paraboloid, thus only the bow shock apex would theoreti- 
cally appear as a spherical central peak brightness distribution, 
where it not for the fact that it will not be visible due to super- 
position with the very infrared bright central star. For the i = 0° 
(side view) the simulation gives indeed A = Ro (i.e. the projected 
and de-projected stand-off distance are the same) and B = Ro V3 
(also here the true apex is observed). 

In Fig. [7] the A and B values, in units of Ro, are given as 
a function of the inclination i. Measuring both A and B would 
then, in principle, directly constrain both Ro and i. Unfortunately, 
measuring both quantities in the detected bow shocks is often 
not possible due to irregularities and incompleteness of the bow 
shock. Also the change of outline with varying i is a small effect, 
thus making the exact inclination angle and stand-off distance 
highly degenerate. For inclination ang les / < 45 ° the cha nge in 
A or B is less than 10% (Fig. [7] and [Mac Low et al.|1991| ). Only 
at large inclinations (or small viewing angles) does the observed 
projected distance A to the bow shock outline differ significantly 
from Ro, the de-projected stand-off distance. Consequently, us- 
ing these relations and the kinematical properties of the mea- 
sured stars and its local ISM, inclination-corrected Ro values can 
be derived by measuring A and/or B together with a calculated 
or assumed i. Thus, first a plausible inclination is derived from 
the space motion, secondly A is measured from the observed im- 
age and thirdly, Rq is obtained from the Ro-A relation in Fig. [7] 
Ideally the star's kinematical properties (proper motion and ra- 
dial motion) would give the heliocentric 3D space motion vector 
of the star, wile the observed bow structure would - indepen- 
dently - yield the heliocentric 3D orientation of the bow, whose 
apex's 3D orientation represents the heliocentric relative space 
motion of the star, with respect to the ISM. Thus, the difference 
of the two would yield the heliocentric 3D ISM flow vector. 
However, in order to circumvent the high degeneracy between 
i and Ro we use the space motion vector inclination angle i to 



constrain Ro from the observed bow shock shape, which thus 
neglects a potential ISM flow (see Sect. |5.1| for a discussion on 
the local LSR ISM velocity). Table [T] gives the de-projected Ro 
obtained via wilkinoid fitting as discussed above for the objects 
which show emission in the shape of an arc or shell, and thus 
potentially trace the dust emission outline of a bow shock. 



4.2. Special cases of "wind-ISM" and "wind-wind" interaction 



In the case of a (nearly) stationary star (with respect to the lo- 
cal ISM and compared to the stellar wind velocity), spherical 
symmetry of the stellar wind-ISM interaction is preserved (e.g. 
Libert et al.|2 007). Nevertheless, there is a relative velocity dif- 
ference between the stellar wind and the ambient medium. The 
interaction of a spherical outflow with external matter leads to 
the formation of a region of compressed material within two 
spherical boundaries ( We aver et al.|1977| ). As before, this region 
consists of a termination shock (where the freely expanding (su- 
personic) stellar outflow is abruptly slowed down by compressed 
material), a contact discontinuity (separating the circumstellar 
and interstellar matter) and the bow shock (external boundary 
at which the external medium is compressed by the expanding 
shell). It is thus a special case of a "wind-ISM" interaction sce- 
nario described in the previous section. The total mass of the 
detached shell would consist of both circumstellar and interstel- 
lar matter. The latter can be estimated from the total hydrogen 
mass that would have been present in a sphere of radius equal to 
the radius of the detached ring (Sect.|3.3|). 



However, the fact that the "external" material of the stel- 
lar wind is interacting with is moving together with the source 
(thus keeping the spherical symmetry) may also suggest that in 
fact the "external" matter is not genuine ISM but rather mate- 
rial remaining from an older mass-loss event, perhaps during 
the RGB phase. The thermal pulse scenario leads to an inter- 
action between a slow and a fast wind, i.e. "wind- wind" inter- 
action. In this case, v* from Eq. [2] is the relative velocity be- 
tween the slow and fast wind. In the rest frame of the star, the 
second wind is computationally equal to a flowing ISM. A de- 
tached dust shell is believed to result from an intense - relatively 
short - episode of enhanced mass-loss rate and wind velocity 
(e.g. initiated by a He- shell flash), whose wind has a higher 
outflow velocity and will thus interact with the previous slower 
wind. The subsequent sharp drop in mass-loss rate and outflow 
velocity for a few thousand years directly following the thermal 
pulse or He - flash will lead to a detached dust shell (e.g.| Qlofsson 
et al.l [T990l [V assiliadi s^Woodl[T993l |Steffen & Schonberner 
2000f |Mattsson et al.|2007j ). The (relative) jump in velocity and 
mass-loss rate from pre-flash to flash-peak is possibly the most 
critical parameter governing the formation of a detached shell 
(Matts son et al.||2U07] ). The latter scenario seems more appro- 
priate for the formation of geometrically thin detached molec- 
ular shells with high expansion velocities showing little inter- 
action with the surrounding medium. Indeed, the detection of 
bow shocks together with a small detached ring gives further 
support to this hypothesis. Another possibility could be that the 
ring represents the (inner) termination shock associated to the 
bow shock located further outwards. Other scenarios, such as 
constant mass-loss in combination with non-isotropic mass-loss 
events and clumpy dust formation have also been detected in the 
form of non-concentric spherical shells (Decin et al. 2011). 



18 



Cox et al.: Bow shocks and detached shells in the far-infrared 



5. Stellar, circumstellar and interstellar properties 

The relation (Eq. [2]) between the stand-off distance of the bow 
shock region and the star's mass-loss rate, M, terminal wind ve- 
locity, v w , peculiar velocity (v*), and the ISM density (nu) is 
powerful in its simplicity. For observed bow shocks around stars 
with known (observed) mass-loss properties and space motion, 
the local ISM density - which is often the most difficult to de- 
termine accurately - can be inferred directly from the measured 
stand-off distance. On the other hand, for stars with known stel- 
lar mass-loss, known space motions and (assumed) ISM densi- 
ties, the stand-off distance can be predicted. In this section we re- 
view the relevant and available stellar properties of the observed 
AGB stars and supergiants as well as the properties of the lo- 
cal ISM. The interplay between these various physical properties 
determines the final shapes and sizes of the wind-ISM interac- 
tion zone, in particular the contact discontinuity as discussed in 
the Sect. |4.1| Adopting the stellar, circumstellar and interstellar 
properties presented in this section for each object in our survey, 
both the predicted ISM densities for observed bow shocks as 
well as predicted stand-off distances for all targets in the survey 
are given in Tables [T] to [3] 

5.1. Stellar distance and (relative) space velocity 

The parallax, proper motion, and the radial local standard of 
rest (LSR) velocities give a direct estimate of the distance, J, 



and the absolute p eculiar space velocities, v* (as per Johns on & 
Soderblo m|1987| ), for the majority of the AGB stars and red su- 



pergiants. The LSR velocities have been corrected for the solar 
motion v n UU, V, W) = (11.1, 12.24, 7.25) kms" 1 ] dSchonrich l 
et al.|201 p] Parallaxes from re-processed Hipparcos data ( |van| 
Leeuwen 2007 ) were taken, when available, if the uncertain- 
ties are less than 20%. For several (distant) stars, the paral- 
lactic distances derived from Hipparcos data have large uncer- 
tainties or have not been obtained at all. For these targets dis- 
tances can be derived, for example, from the (pulsation) period- 
luminosity (K) relation, and observed apparent K magnitudes. A 
generic Galactic P-K relation valid for both oxyge n- and carbon- 
rich (Mira) v ariables was established by e.g. Whitel ock et al. 
([20001 [2008] ): M K = -3.51 [log P - 2.38] - 7.25. |Vlemmmgs 
|et al . | ( [2003] ) obtained VLBI parallaxes for four Miras that were 
more consistent with distances derived via the P-K relation than 
those obtained with Hipparcos. Tables [T]to[3]give the derived pe- 
culiar velocities and peculiar absolute proper motions with their 
associated position angles. 

Radial (LSR) velocities are taken from CO line surveys 
where available. Prim ary sources are pe Beck et al. ([2010 ) and 
Men zies et al. ( 2006| ), which have good agreement for sources in 



common to both studies (typically the differences are less than 
3 kms" 1 ). 

At this point we assume that the relative peculiar velocity be- 
tween the ISM and the star, v*, is determined entirely by the star 
space velocity with respect to the local standard of rest (LSR) 
(i.e. the stationary ISM). In other words, we assume that there 
is no flow of the ISM itself. This could be a simplification for 
some cases where the ISM may have an appreciable flow veloc- 



ity as matter is being blown away by super-bubbles (see |Ueta 
et al. 2009a b for the case of a Ori). 

To test this, we estimate the loc al LSR ISM velocity 
(v r ,v/,Vfc) from the galactic rotation (Oorf| [T928] |Feast & 



1 Note that other common values adopted for the solar motion from 
Mihalas & Routly ('1968'), Dehnen & Binney (1998 ), or Reid et al 

J. „ I -: * * . ' . .„ U . ' \' = 



2009 ) give results for v* that differ slightly, by a few km s 1 . 



Whiteloc k|1997| ) at the position of the stars and correct v* to ob- 
tain the local v*_ism- For the majority of the objects the v*_ism 
is within ±20% of v*, thus decreasing/increasing the predicted 
stand-off distance by the same factor. Only in a few cases this 
correction gives a significant change in peculiar velocity; for 
U Cam, V Hya, T Lyr, and [i Cep this gives relative space veloc- 
ities that are approximately a factor of 1.5 higher, and for YPav 
a factor of 2 lower. 

The average space velocity, v*, for our sample of stars is 
-36 kms" 1 (Tables to |3}, consistent with the average velocity 
of -30 km s" 1 for Galactic AGB stars ( [Feast & Whitelock|2000[ ). 



5.2. Mass-loss rate 

The mass-loss rate of evolved stars is important for the enrich- 
ment of the ISM. It also directly affects the shape and size of the 
stellar wind-ISM interaction (Sects. [4] and [6]). Mass-loss rates 
can be estimated directly via modeling of observ ed CO line pro- 
files (e.g. Knapp et al. 1998[ Groenewegen & de Jong |1998| 



Groenewegen et al.|2002[|De Beck et al.|2010j ). In general the de- 
rived (gas) mass-loss rates depend on the line strength, terminal 
gas expansion velocity (v^), and distance (d 2 ). The average gas 
mass -loss rate from -300 Galac tic carbon stars is 1.1 x 10" 5 M© 
yr" 1 (Groenewegen et aL|2002|). The latter should probably be 



reduced by an order of magnitude since new Hipparcos results 
show that the distances adopted by Groenewe gen et al.] p002| ) 
are over-estimated by a factor of two to four. Using CO-derived 
mass-loss rates, an empirical relation has been established be- 
tween the mass-loss rate and luminosity variability period (e.g. 
|De Beck et aT]|2010| ): log(M) = -7.37 + 3.42 x 10" 3 P (for P 
< 800 days) and log(M) = -4.46 (for P > 850 days), with M in 
units of M n yr -1 . However, the scatter on this relation is rather 



large with a typical uncertainty of a factor of 10 (see De Beck 
|et al.||2010] for det ails). Similar results are obtained for carbon 
Mira variables by Groenewegen & de Jong ( 1998| : log(M) = 
4.08 log P - 16.54. Adopted mass-loss rates - based on CO ob- 
servations or luminosity period - are included in Tables [T] to [3] If 
neither M nor period are available, we adopt a generic value of 
5 x 10" 7 M yr -1 . Uncertainties in M arise predominantly from 
inaccurate distance determinations for distant stars. 



5.3. Stellar wind velocities 

For the (terminal) wind velocit y, v W9 we adopt the ter minal ve- 
locity of the CO envelope ( e.g. |De Beck et al.||20T0t Tables [T] 
to [3] and references therein). Groeneweg en et al. ( |2002| ) find an 
average expansion velocity of 18.7+6.1 kms" 1 for a sample of 
330 carbon stars. For a set of 24 oxygen-rich and 13 carbon-rich 
AGB stars, |De Beck et al.| ( |2010| ) find mean v w values of 14.5 and 
15.4 kms , respectively. These values are 4-10 kms -1 higher 
than those of |Ramstedt et al . (2006) for a set of 77 oxygen- 
rich and 61 carbon-rich stars. For carbon-rich stars, this average 
3 kms -1 lower than the one found by |Groenewegen et al. 



is 



( 20021). If no CO terminal velocities are available, we adopt a 
generic value of v w = 15 kms -1 , for both oxygen- and carbon- 
rich stars. 



5.4. Local interstellar medium density 

The ambient (uniform) mass density of the ISM is defined as 
Pism = y^H ^h, with //H = 1.4, the mean nucleus number per 
hydrogen atom for the local medium, m H is the mass of the hy- 
drogen nucleus, and nu the interstellar hydrogen nucleus density. 
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Although the ISM is far from uniform, there is a general depen- 
dence between the interstellar matter (space-averaged) volume 
density, ttism, and z, the distance from the Galactic plan^] 

The different phases of the ISM have different average den- 
sities, scale-heights and filling factors (see e.g. |Savagej|1995| 
Boulanger|2Q0T]|Ferriere|2001t|Wooden et al.|2004[|Cox|20051 ). 
Dickey & Lockman (1990) used radio observations of Hi to de- 
rive the structure of the CNM, WNM and WIM (warm ionised 
medium):. 



n H (z) = 0.39 e 



( 127pc ) 



+ 0.11 e v3i8 P c) + 0.06 e ^ 



(4) 



Bas ed on UV observations o f Hi and H 2 , Sav age et al. ( 1911) 
and piplas & Sava"ge ( 199 A) arrive at a similar result for the 
WNM as|Dickey & Lockman] ([1990]). An alternative relation is 



given by |Loup et al 



( 1993| ) w hich is based on Spitzer ( T978J) 
2cm~ J ) and |Mihalas & Binney||l98l| (scale 



(density ?7h(z=0) = 
height of 100 pc): 

n H (z) = 2fie~^ (5) 

These two relations are shown further below in Fig. pT] (Sect. [7]). 
Thus, we obtain some first insights into the local ISM density 
for our objects (Tables [T] to [3]). Evidently, one should be cau- 
tious using these relations for specific cases as there is struc- 
ture in the ISM on all spatial scales, with different phases that 
have different filling factors. For example, diffuse to molecular 
clouds (hh ~ 10 - 1000 cm -3 ) have an order of magnitude lower 
filling factor than the low density WNM and WIM surround- 
ing it. Thus, statistically about 10% of the objects in our survey 
could be moving through a denser medium than inferred from 
the above relation for the WNM. Generally, these relations do in- 
dicate that Tin drops as one moves away from the Galactic plane, 
leading to - on average - lower volume densities of the ISM. We 
use Eq.[5]to get a first estimate of the local ISM density for the 
objects in this survey (Tables [T] to [3]). Accurate measurements of 
the local ISM density for all stars in our sample would be desir- 
able, but are currently - to the best of our knowledge - unavail- 
able. In fact, we will show that for certain cases the observations 
of bow shocks can potentially be used, via the measurements 
of the stand-off distances, to derive estimates of the local ISM 
density. 

5.5. Circumstellar chemistry, spectral type and binarity 

In addition to the above stellar properties that have a direct im- 
pact on the theoretical stand-off distance via Eq. [2] other pos- 
sibly pertinent information such as spectral types of the central 
star and dominant circumstellar chemistry (oxygen versus car- 
bon rich), and binarity is also included in Tables [T] to [3] If and 
how these stellar attributes could affect the occurrence and shap- 
ing of bow shocks is discussed further in Sect. [7] 

6. Hydrodynamical models of interaction between 
the slow stellar wind of a moving star and the ISM 

Hydrodynamical simulations offer the opportunity to explore 
the effect of varying physical properties of either the stars (e.g. 



2 z can be expressed as z(pc) = d sin b + z Q , with d the distance, b 
the Galactic latitude, and z Q the sun's vertical displacement from the 
Galactic plane. The exact value of z Q depends strongly on the assumed 
underlying Galactic model and/or the observational data selection cri- 
teria. However, the v alue is converging towards - 15 - 20 pc ( 13+7; 
Brand & Blitz 199 3; |Humphreys & Lars"enp995[ 20.5+3.5 pc; |Reed| 
2006 19.6+2.1 pc; |Joshi|2007| 17+3 pc). We adopt a value of 15 pc. 



mass-loss rate, wind and (relative) space velocity) and/or the 
ISM (density, temperature) in a coherent systematic way. For 
example, |Steffen & Schonberner| ( [2Q00| ) and |Libert etaL]pOQ7| ) 
used hydrodynamical simulations to show that a brief episode 
of increased mass-loss rate could give rise to an expanding, ge- 
ometrically thin shell. |Wareing et al.| ( |2006b| ) applied numeri- 



cal simulations of a two-wind model t o explain the observed 
structure around RHya. Simulations by Warei ng et al.| |2QQ7a) 
and |Wareing et al.| ( |2007b| ) indicate that a higher mass-loss rate 
(for similar ISM density and space velocity) will result in more 
pronounced Kelvin-Helmholtz instabilities (see their Fig. 1). 
Wareing et al.| (|2006a) obtain bullet-sha ped emission structures 
simulating a PN (their Figs. 3 and 5). Villav er et al. (2003 ) in- 
clude time-dependent mass-loss in simulations of wind-ISM in- 
teraction, leading to time-dependent stand-off distances. 

One particular strength of hydrodynamical simulations is 
that they allow to study the formation, growth and dissipation 
of fluid instabilities - like Rayleigh-Taylor (RT) and Kelvin- 
Helmholtz (KH) instabilities - important in wind-ISM shock in- 
teractions. RT instabilities occur when a dense, heavy fluid is 
accelerated by a light fluid. Normally two plane-parallel fluid 
layers are meta- stable, but a small perturbation can destroy this 
delicate equilibrium. This manifests itself as so-called inter- 
penetrating "RT- fingers". Such instabilities could be quenched 
by a restoring force such as a magnetic field, th us preventing 
these instabilities to grow (Chandrasekhar 1961). The KH in- 
stability results from a velocity shear between two fluid lay- 
ers. This flow of one fluid over another will induce a centrifu- 
gal force which leads to changes in pressure which amplifies 
the ripple. Together with a RT instability, the KH instability will 
form structures in the shape of mushroom caps on the end of the 
RT fingers. The KH time-dependent turbulent eddies can form a 
complex structure arising in a steady flow, if they become large 
enough to influence that large-scale morphology of the shocked 
gas. Other i nstabilities that can occur are the non-linear thin shell 
instability ( Vishniac 1994) and the transverse acceleration insta- 
bility ([Dgani et al.|1996|) as shown numerically b y e.g. |Blondin 



& Koerwer| ( |1998| ) and |Comeron & Kaper| ( |T998l ). 



6.1. Hydrodynamical simulations 

Here we present a series of seven simulations of the interac- 
tion between the ISM and the circumstellar medium of moving, 
evolved stars in order to find out how the morphology of the 
bow shock varies with the various stellar wind and ISM param- 
eters. The different pertinent parameters for the simulations are 
summarised in Table [6] For our hydrodynamical simulations we 
use the MPI-AMRVAC code ( [Keppens et aLl[20TT] ). This code 
solves the conservation equations of hydrodynamics on an adap- 
tive mesh (AMR) grid. We use a 2-D cylindrical grid in the r-z 
plane. The basic resolution is set at 80 grid points per parsec, but 
allows four additional levels of refinement, each doubling the ef- 
fective resolution. This gives us a maximum effective resolution 
of 1280 grid points per parsec. Since some of the models require 
a larger physical space we use grids of different sizes for those 
simulations. In those cases we increase the number of grid points 
on the basic level to maintain the same resolution. The expand- 
ing stellar wind is inserted by filling a small circle with wind 
material. The motion of the star is handled by giving the ambient 
medium around the star a velocity in the z-direction. Therefore, 
we are simulating the wind interaction in the frame-of-reference 
of the star. Since all models are 2-D and the simulated space lies 
along the direction of motion of the star, the snapshots (Figs. [8] 
through [TO]) show projections for i = 0. 
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Table 6. Parameters for the moving star simulations. 



label 


description 


M 


M dust /M gas 




^ISM 


^ISM 


physical space 


basic grid 


R a 


d b 






(M© yr 1 ) 


(%) 


(km s v ) 


(cm 3 ) 


(K) 


(pc) 




(pc) 


(pc) 


A 


basic modeF 


10" 6 


1 


25 


2 




1.5x1 


120 x 80 


0.26 + 0.01 


0.35 


B 


high M 


10- 5 


1 


25 


2 




2x2 


160 x 160 


0.59 + 0.03 


0.73 


C 


low M 


10- 7 


1 


25 


2 




2x1 


160x80 


0.090 + 0.001 


0.10 


D 


low dust 


10" 6 


0.1 


25 


2 




2x1 


160x80 


0.33 + 0.03 


0.33 


E 


high v* 


10" 6 


1 


75 


2 




1.5x1 


120 x 80 


0.11 +0.02 


0.13 


F 


low nisu 


10" 6 


1 


25 


0.2 




3x2 


240 x 160 


0.76 + 0.08 


0.92 


G 


warm ISM 


10" 6 


1 


25 


2 


8 000 


2x2 


160 x 80 


0.31 +0.06 


0.45 



Notes. (a) Estimated. Large-scale instabilities make R (contact discontinuity) a time-dependent property. In particular, for simulations E and G 
the mixing is very efficient, eliminating in effect the contact discontinuity altogether. R s , the location of the forward shock of the shocked gas 



region, where the ISM transition from unshocked to shocked gas can be more accurately measured 



(c) - 



15 km s for all models. 



To the basic conservation equations we have added the effect 
of optically thin radiative cooling. This is necessary since some 
of the shocks are strongly radiative, which changes the morphol- 
ogy of the shocked gas. In the case of simulation G, where the 
ISM has a high temperature, we put a lower limit on the cool- 
ing in the ISM, so that in the ISM the temperature does not fall 
below 8 000 K. More importantly, we have also added a (sim- 
plified) dust component (van Ma rie et al.||2011b| ). This is done 
by using a two-fluid approximation, with the dust represented as 
a gas without internal pressure. For simplicity we only include 
dust in the wind and neglect the dust component of the ISM. This 
gives us a set of five partial differential equations (not counting 
the vector components) for the pressure balance for the gas and 
dust. For the gas we have conservation of mass, momentum and 
energy, while for the dust we only have conservation of mass and 



momentum. The appropriate equations are given in |van M arie 
et al.| ( |20lTb| ). Radiative losses depend on the hydrogen and elec 



tron particle densities (derived from p assuming full ionization 
with hydrogen mass and involve a temperature dependent cool- 
ing curve A(T ). This coolin g curve has been calculated with the 
CLOUDY code (Ferla ndet al.|1 998) and includes radiative losses 
through IR radiation from the dust. We also include a drag force 
linking the gas and dust, which is derived from a combination of 
Epstein's drag law for the s ubsonic r egime and Stokes' drag law 
for the supersonic regime ( |Kwok|1975] see also | van Marie et alT 



201 lb). The drag force depends on the dust particle density and 
radius (0.005 /mi), the velocity differe nce between g as and dust, 
as well as the thermal speed of the gas ( Kwok 1975 ). For further 
details we refer the reader to | van Marie et al.| ( |201 lb\ . 

Since our equations are purely hydrodynamical, we don't 
take into account magnetic fields or the effect of radiation on the 
dust particles. We also neglect destruction and creation of dust 
particles. Radiative processes are not included and thus the dust 
temperature cannot be treated appropriately (currently only col- 
lisional heating of dust could be accounted for). Also the gas 
is purely heated by collisions, and cooled radiatively; photo- 
ionisation is not included. Introducing radiative cooling by dust 
will cause the bow shock to become thinner and more unstable. 

Adding destruction processes such as collisional heating or 
UV irradiation could, in principle, be included, but at a high cost 
in the required computing time. Both mechanisms are not ex- 
pected to be very efficient as grains are not easily destroyed, 
because they are very effective radiators (emitting all heat im- 
mediately in the infrared) and the shocks discussed here are not 
very strong, and the radiation field would not be strong enough to 
destroy grains. Furthermore, dust can also be created in the high 
density shock regions, thus lowering t he effective dust destr uc- 
tion rate. Further details are given in van Mari e et al. ( 201 lb| ). 




Fig. 8. Top panel (A): Gas density in g/cm 3 (right) and dust 
grain particle density in cm -3 (left) for simulation A after 
1.37xl0 5 years; basic model. The bow shock is smooth, but 
at the contact discontinuity both Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities are visible. Bottom panel (B): Similar to 
top panel, but for simulation B after 5.0xl0 4 years; high M. The 
instabilities are now primarily of the Rayleigh-Taylor type. Due 
to the stronger wind the bow shock lies further from the star. 



The input and grid parameters are summarised in Table[6] We 
start with a basic model (simulation A), whic h has input param - 
eters based on the observations of a-Orionis ( |Ueta et al.|2008j ). 
Using this model as a starting point, we vary individual parame- 
ters to investigate the effect on the morphology of the bow shock. 
Note that for a stationary star, v* = km s" 1 , the shock between 
the wind and ISM will drive a spherical shell of wind material 
sweeping up the ISM material, thus leading to detached shells / 
rings (not shown). 
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E - fast moving star 



Fig. 9. Similar to Fig. [ij but for simulation C after 
1.0x10 s years(top) and simulation D after 5.0xl0 4 years (bot- 
tom). Top panel (C): Because of the weaker wind the bow shock 
is very close to the star. The bow shock morphology is com- 
pletely stable. The shocked gas fills up the cavity behind the star 
due to low ram pressure of the wind, and thermal pressure of 
the shocked gas, because the shocks are almost completely adi- 
abatic. Bottom panel (D): The bow shock is more conical than 
for the simulations with a stronger dust component and shows 
more instabilities. The instabilities are primarily of the Kelvin- 
Helmholtz type. 



6.2. The morphology of the bow shock 



As shown in Figs. [8] through 10 the collision between the ex- 
panding wind and the moving (in the rest-frame of the star) 
ISM creates a bow shock. The location and morphology of the 
shocked gas depend on the exact input parameters. These figures 
show the gas density and dust particle density in the circumstel- 
lar environment after the bow shock has reached its (semi-) per- 
manent equilibrium distance from the star. Typically, this takes 
between 50000 and 150 000 years, depending on the parame- 
ters of the simulation. Our basic model (Simulation A, Fig. [8} 
shows a standard bow shock morphology: a forward shock sep- 
arating the unshocked and shocked ISM, a wind termination 
shock that separates the free- streaming wind from the shocked 
wind, and between them a contact discontinuity separating the 
shocked wind from the shocked ISM, with the temperatures in 
the shocked wind typically lower than in the shocked ISM due 
to the higher density of the former. This pattern is repeated in 
all subsequent models, with the exception of the fast moving 
star (Simulation E, Fig.[l0]) and the high ISM temperature cases 
(simulation G, Fig.[l0]), both of which show a chaotic bow shock 
with strong turbulent behaviour. 

Behind the star a relatively structureless - empty - region re- 
mains. This is only filled by the free- streaming wind of the star. 





1 pc 



Fig. 10. Similar to Fig. [8j but for simulation E after 
6.5xl0 4 years (top), F after 5.0xl0 4 years (middle), and G after 
1. Ox 10 5 years (bottom). Top panel (E): The entire shocked gas 
region is completely unstable due to a combination of radiative 
shocks and the high ram-pressure from the stellar motion in the 
ISM. The high ISM ram pressure brings the bow shock close to 
the star and the region behind the star is filled up with shocked 
gas due to the turbulent motion of the bow shock. Middle panel 
(F): The low density ISM results in a very extended bow shock 
region. Bottom panel (G): Due to the high ISM temperature, the 
radiative cooling in the shocked ISM is relatively ineffective. 
Since the shocked ISM maintains a high thermal pressure, due 
to the enforced minimum temperature of 8 000 K, the shocked 
gas fills the void behind the moving star. 



Note that in our simulation the wind velocity is always less than 
the stellar velocity. Exceptions occur if the shocked gas region 
has a high thermal pressure, which pushes it into the empty re- 
gion (Simulation C, Fig. [9] and simulation G, Fig. 10), or if the 
shocked gas is very turbulent (Simulation E, Fig.|10|). 

The location of the shock can be approximated by the wind 
and ISM conditions as described in Eq. [2] In Table [6] we give 
the location of the contact discontinuity (Rq) as well as the loca- 
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tion of the forward shock of the shocked-gas region (R s ), where 
the ISM transition from unshocked to shocked gas is located. 
In some cases Rq has to be approximated since instabilities can 
make it a time-dependent value. For simulation A (the basic 
model), the analytical approximation puts the contact disconti- 
nuity at 0.42 pc. In our simulation the position of the contact dis- 
continuity deviates from this prediction, lying at about 0.26 pc 
(while the forward shock lies further ahead at -0.35 pc). The 
difference can be explained as a result of radiative energy loss 
in the shock. The scaling, as described by Eq. [2] appears to be 
quite accurate. For example, for an order of magnitude decrease 
in ISM density (simulation F, Fig. [T0|) or increase/decrease of 
mass-loss rate (simulation B, Fig. |8f7 the location of the bow 
shock moves away/closer with a factor of about VlO. 

The dust grains, which start with the same velocity as the 
wind, are initially carried along. When the wind reaches the 
termination shock, the gas slows down abruptly. Because they 
are not subject to the shock, the dust grains keep moving, but 
are now subject to drag force due to the difference in velocity 
with the gas. This causes them to slow down over time. As a re- 
sult they tend to pile up at the contact discontinuity ( [van M arie 
et al.||20lTb| ). Once the grains reach the contact discontinuity, 
they tend to follow the local instabilities and are eventually car- 
ried downstream. Those grains, which crossed the contact dis- 
continuity, will usually go downstream faster, since in most of 
our simulations the shocked ISM has a higher velocity than the 
shocked wind. The influence of the grain size on the gas-dust 
interaction is investigated in van Mar ie et aL] ( |201 lb| ). 



abatic, the shocked-gas region is wide and its thermal pressure 
pushes material into the area behind the star, rather than leaving 
it empty. 



6.3. Instabilities in the bow shock 

In most of our simulations, instabilities occur along the contact 
discontinuity. These instabilities consist of a combination of RT 
and KH effects (see above). In our case, the (relatively) low- 
density- shocked ISM exerts a force on the denser shocked wind, 
which leads to the formation of RT fingers. As a result, the wind 
material starts to flow into the shocked ISM. Here the wind ma- 
terial is subject to a sheer-force due to the relative motion of the 
star with respect to the ISM. These instabilities are enhanced by 
the dust (albeit this is a secondary effect), which tends to con- 
tinue in a straight line in an attempt to cross from the shocked 
wind into the shocked ISM, dragging the gas along with it. This 
is best observed in the case of simulation B (Fig. [8]). The ve- 
locity difference between the shocked gas layers causes KH in- 
stabilities to develop wherever the RT instabilities cause a local 
displacement of gas. These lead to the characteristic "cyclonic" 
features (see Figs.[8][9]and[T0]). The presence of dust slows down 
the formation of KH instabilities, since the circular motion of the 
KH instabilities has difficulty overcoming the inertia of the dust 
grains. Therefore, they are most visible in simulation D (Fig. [9]). 
Also, the KH instabilities can form more easily if the low density 
region (the shocked ISM) is extended, so that the instabilities do 
not hit the forward shock (simulation G, Fig. 



10) 



In the case of the star with a low mass-loss rate (simula- 
tion C, Fig. [9]), such instabilities are absent. The low mass-loss 
rate leads to a reduced density in the shocked wind. As a re- 
sult, the density difference between the shocked wind and the 
shocked ISM is small, which reduces the RT effect. Also, due to 
the low density, the radiative cooling, which scales with th e den- 
sity squared, is less effective ( |van Marie & Keppens|201 1 ). This 
is important because radiative cooling tends to favor the forma- 
tion of small, high density clumps, which in turn can serve as 
a start of other instabilities. Because the shocks are nearly adi- 



The fast moving star (simulation E, Fig. 10) has a far more 
irregularly shaped bow shock. In this case the shocks on both 
sides are highly radiative, leading to a very compressed shocked 
gas region. This, combined with the strong ram-pressure from 
the ISM, causes the entire shell to become unstable, leading to a 
ragged form which will change considerable over time as local 
instabilities grow to a size where they dominate the entire struc- 
ture of the shell. This turbulent motion also causes the gas (and 
the dust grains) to end up behind the star. 



7. Discussion 

7.1. Comparison between observed and hydrodynamical 
bow shock morphology 

The simulations presented in Sect. [6] only cover part of the pa- 
rameter space spanned by all stellar, circumstellar and interstel- 
lar properties. In particular, the simulations focus on the effects 
of increasing the star's peculiar velocity with respect to the ISM, 
changing mass-loss rate, increasing the ISM temperature and 
lowering the ISM density. 

One evident issue revealed by the simulation is that the 
shocked gas region can have considerable spatial extent (this oc- 
curs when the shocks are adiabatic - as in most cases presented 
here - as opposed to radiative). In the approximation of efficient 
cooling, the shock region will be thin and unresolved; however, 
the simulations show that this is not necessarily valid for all 
cases and can lead - depending also on the exact cooling law - 
to extended interaction zones. Depending on the exact morphol- 
ogy and the temperatures in the shocked gas region, one could 
- theoretically - be observing either the forward shock, reverse 
shock or contact discontinuity. However, the current far-infrared 
observations are most sensitive to cold dust grains. Due to the 
limited spatial resolution, it is not clear which region is actually 
represented by the far-infrared emission. It may represent the 
contact discontinuity (i.e. the entire unre solved shocked r egion), 
the bow shock, or the termination shock ( Libert et al. 2007 ). Our 
simulations in Sect. [6] show that most dust grains pile-up at the 
contact discontinuity where they tend to follow the local instabil- 
ities which carry the grains downstream. Thus, the simulations 
suggest that the observed far-infrared dust emission primarily 
traces the contact discontinuity. On the other hand, if grains are 
destroyed in the shock they will be found in the unshocked stel- 
lar wind region and thus far-infrared dust emission would delin- 
eate the termination shock (and the contact discontinuity if the 
shock region is physically thin enough). This might prove useful 
in explaining the spatial offset between observed UV emission 
and far-infrared emission in the bow shock of CWLeo (Ladjal 
|et al.||2010| ). Further work on the destruction and formation of 
dust grains (c.q. alteration of the dust size distribution) in shocks 
is warranted. As stipulated above, the simulations do not include 
radiative transfer and can thus not provide appropriate dust tem- 
peratures needed to simulate the infrared emission. 



RT and KH instabilities 

The presented PACS infrared observations reveal RT and KH 
instabilities in astrophysical bow shocks. RT fingers can be seen 
in the shock region of R Scl. KH "wiggles" or density knots are 
more common, and can be seen in the bow shock regions of, for 
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Table 7. Comparison between observed stand-off distances de- 
rived from Herschel/PACS maps and previous studies with 
IRAS, Spitzer & AKARI. The predicted values (this work) are 
given as a function of hh (cm -3 ). The local interstellar ISM den- 
sity, hh, is estimated from the predicted Rq/ ^JKh and the ob- 
served Rq and compared to n H given by Eq.[5] 



Object 




De-projected R 




Ref 




n H 
[cnT 3 ) 




this work 


literature 


observed (') predicted (') 


(') 


Eq. 


2\ Eq.[5] 


oCet a 


1.2 


0.7/ V"// 


~3 


(1) 


0.4 


1.1 


a Oh 


5.0 


10.6/ 


~4 


(2) 


4.6 


1.9 


CW Leo 


6.6 


3.9/ 


5.9 


(3) 


0.3 


0.7 


RHya 


1.6 


2.4/ 


1.6 ±0.1 


(4) 


2.5 


0.8 


RCas 


1.5 


2.2/ V^" 


1.4 ±0.1 


(5) 


2.1 


1.8 



Refe rences. References: (1) Raga & Canto 20 08) (2) 
Ueta et al. 200 8] Decin et a l. i n preparation; (3 ) |Ladjal e t al. 2010 



|Sahai & Chronopoulos|2010l (4) |Ueta et al.| ( [2006] ); (5) IUeta et al.l20T0l 

Notes. (a) Note that for o Cet the observed features delineate most likely 
the termination shock and not contact discontinuity, this measure is thus 
a lower limit to R , which will be larger than \.2'(e.g. more like 2' as 
observed in the UV) leading to a lower density of n H = 0.2 cm -3 . 



example, UUAur, RHya, XPav, EPAqr, //Cep, RLeo, RT Vir, 
XHer, andV1943Sgr. 

Almost spherical bow shocks spanning a large azimuthal an- 
gular range correspond closest to simulations A, B and F, the 
nominal, high M and low cases, respectively. For low 
(simulation F) the KH instabilities are more pronounced, while 
the RT instabilities are more evident in the high M case. The 
smooth, spherical bow shocks of a Ori, UU Aur, and X Pav are 
similar to the morphology of the simulation A (nominal case; 
Fig. [5]). UUAur and XPav show back flow emission also seen 
in the simulation. The "wiggles" in the bow shocks of R Hya 
and RLeo correspond to those seen in simulation F (Fig. [10| ) 
which contradicts the high ISM density derived for RLeo (see 
Sect. 7.5). Both XHer and V1943Sgr have high space veloci- 
ties leading to morphologies much like that seen in simulation E 
(Fig. [TO}. 

Further comparisons are hampered on the one hand by the 
(unknown) projection effects in the observations and on the other 
hand by the, currently, limited parameter space covered with the 
few simulations at hand. A more extensive, multi-parameter grid 
of simulations and/or a detailed modeling of individual cases 
will be required to match observations with simulations to con- 
firm whether or not the latter accurately predict the former. 



emission. The proper motion for Class III objects is not consis- 
tently aligned with, for example, shell features. In a few cases 
(TT Cyg, UAnt) there is some stronger shell emission in the 
direction of the proper motion, but this is not observed for the 
other detected rings. 

There is no clear indication for a dependency on the inclina- 
tion angle, i. All five classes include sources with a range of in- 
clinations between almost 0° up to 60-70° . Except for VY UMa 
(/=-l°), there are no objects with small inclination angles among 
the "eyes" (all have i > 28°), whereas the " fermata' ', "rings" 
and "irregular" all have a more uniform distribution over incli- 
nation angles. However, we stress that at this point the sample is 
too small to derive any further conclusions. 

For about one third of the "non-detection" class, the proper 
motions are not known and thus their space velocities are, at 
best, lower limits based on radial velocities. This is in part a 
distance bias. The objects with non-detections are on average 
further away and consequently their proper motions have been 
too small to be measured. This means also that the predicted 
stand-off distances are upper limits. 



7.3. Binary interaction, circumstellar chemistry, and magnetic 
fields? 

Binary interaction 

There is little evidence for binary interaction in our sample. Even 
though for several stars in our sample, there is solid evidence 
for close companions (see Tables [T] to [3}, none show suspicious 
features in their bow shock interfaces. However, many of these 
appear in the category "irregular", which could indicate that 
such wind-wind-ISM interaction is only visible on small spatial 
scales. It is difficult to arrive at any firm conclusion as the num- 
ber of binaries in the survey is low. Excluding the optical "vi- 
sual" binaries, as these are not yet confirmed, out of the 25 ob- 
jects in the "fermata" class only one is a known binary (Mira) 
which indeed reveals a peculiar shape of the interaction zone 
(Mayer et al. 2011). Three others, W Aql, EP Aqr and #Aps are 
solid candidates for binarity. None of the "rings" are binaries. 
Noteworthy, 3 out of 6 of the "irregulars" are confirmed bi- 
naries, and 2 of the 6 "eyes" are potentially binaries. Finally, 
also 6 out of 30 non-detections are binary systems. Including vi- 
sual binaries adds 3 binaries to each of the "fermata", "eyes", 
and "non-detection" classes, and 4 binaries to the "rings". In 
this case, 5 out of 7 "eyes " show evidence for binarity, whether 
this means there is a connection between binarity and the "eyes " 
morphology can not be confirmed nor excluded at this point. 



7.2. Proper motion and inclination 

In most cases the proper motion vector (if known) points roughly 
(within -20° -30°) in the direction of the observed (Class I) bow 
shock apex. Exceptions are \i Cep and o Cet. One possible ex- 
planation is that, like o Cet, \i Cep is also a long-period (tens of 
years) binary, for which the present-day proper motion is not 
pointing in the same direction as it was when the bow-shock 
matter was expelled. This conjecture needs to be further tested 
by further observation and modeling of these systems and their 
observed bow shocks. 

For the majority of the Class II objects, the proper motion 
direction is aligned roughly along the symmetry axis (perpen- 
dicular to the arc-star-arc line). Only for VPav the proper mo- 
tion points to the north-west arc. For Class IV the proper mo- 
tion direction coincides with the location of extended irregular 



Oxygen-rich versus carbon-rich chemistry 

One key property of the AGB stars and red supergiants in our 
sample is their circumstellar envelope chemistry. This can be ei- 
ther carbon-rich (the dust is predominantly made up of amor- 
phous carbon, graphite and silicon carbide) or oxygen-rich (the 
dust is composed of silicates and oxides). Although we do not 
find a strong dependence on the presence versus absence of bow 
shocks and rings with respect to chemistry, there does appear 
to be a distinction between the shape of the extended emission 
found around carbon and oxygen-rich objects. The "fermata" 
and "irregular" classes include a high fraction, 22 out of 31, 
of O-rich stars. The "eyes" and "rings" , on the other hand, in- 
clude only few, 3 out of 22, O-rich stars. It seems that at least 
the "ring" structure is typical of C -rich chemistry, as already 
pointed out by Olofs son et al. ( 2000| ). A further argument favor- 
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ing such a link between "ring" structure and C-richness is that 
all "fermata" stars also displaying a "rings" structure are C- 
rich stars (they are noted with an asterisk in Table [T] and are also 
listed in Table [5]). 

Could the circumstellar envelope structure be related to 
chemistry through differences in mass loss rates? It has been 
shown ( Knapp 1985) that oxygen-rich stars have a wider range 
of mass-loss rates than C-rich stars for which the M is globally 
more homogeneous. On average, mas s loss rates of C s tars are 
higher than those of MS and S stars ([Guandalini | |20 1 0| ) . Apart 
from the fact that the mass-loss rate standard deviation of O-rich 
stars is indeed found to be somewhat larger than that of C-rich 
stars, no clear dichotomy emerges between the mass loss rates 
of O-rich stars and C-rich stars in the present sample. 

Furthermore, the higher emissivity (Draine & Lee|1984 Li 



|& Draine|2001| ) and lower gas-to-dust ratio of 160 for oxygen 
rich dust co mpared to the higher ratio of 40 for carbon-rich 
dust (see e.g. |Knapp|1985l|Heras & Hony|2005) will, for O-rich 
dust, give rise to brighter infrared emission for a given total dust 
and gas mass. One could then argue that the O-rich stars mass 
loss is more easily detected even with a patchy structure. This 
could perhaps explain why O-rich stars are detected preferen- 
tially among both "irregulars" and "fermata" classes. 

An alternative explanation is that the "Ring " (and possible 
the "Eyes ") morphology, are explained through (i) the interac- 
tion between a fast wind sweeping out matter from a previous 
slower wind, or (ii) a phase of drastically enhanced mass loss, 
e.g. caused by a thermal pulse. If this explanation holds, and if 
one assumes that most thermal pulses events are followed by a 
third dredge-up (TDUP) episode, then the discriminant charac- 
ter would not be O-rich or C-rich, but instead pre-TDUP or post- 
TDUP (i.e., thermally-pulsing (TP) AGB) stars). All intrinsic S- 
type stars, C-type stars and technetium-rich M stars in Tables [T] 
and [3] are TP- AGB stars. In fact, all "Ring" stars are TP- AGB 
stars. 



Stellar and interstellar magnetic fields 

Some wind-ISM models include the effects of an interstel- 
lar magnetic field on the shaping of the bow shock region. 
Heiligman (1980) predicted that the ISM magnetic pressure 
could deform the outer shells of very large spherical halos of 
PNe into "lemon" shapes, with the axis of symmetry inclined 
to the ordered ISM magnetic field direction. |Soker &T) gani 
(1997]) and |DganT| ( |1998| ) elaborated on this scenario predict- 
ing different bow shock shapes depending on the relative magni- 
tudes of the star's space velocity, v*, the wind velocity, v w , and 
the Alven speed, va- The double arc objects (Class II) such as 
VY UMa and AQ Sgr could perhaps provide evidence for this 
scenario, although both are at relatively high latitudes and their 
symmetry axes are not aligned with the Galactic plane (assum- 
ing that at these locations the magnetic field is also parallel 
to the Galactic plane). It is also possible that the star has its 
own magnetic field and thus may exhibit axis-symmetric mass 
loss. Preliminary simulations (private communication A. -J. Van 
Marie) indicate that double arcs as well as jets (tentatively iden- 
tified for VY UMa and W Hya) could be formed by a (rotating) 
star with a magnetic field. 



or assumed expansion velocities, distances and angular sizes. 
Ranging from the youngest shells like U Cam or R Scl with ages 
of the order of 1000 years and much older ones like AQAnd, 

Y CVn or UX Dra with ages of a few 10 000 years we see also a 
significant trend in the detectability in fossil mm-CO. Whereas 
all objects with shell ages below some 1000 years are prominent 
in CO, no detached shell older than 10000 years was success- 
fully detected yet. This could be understood in terms of photo- 
dissociation by the interstellar radiation fie ld in the older ex- 
tended shells with sizes of 0.2 pc or more ([Kerschbaum et al. 
|20T0| ). |Arimatsu efaTI ( [20lT| ) find that the detached shell of 
U Ant detected with PACS is not the outermost one found with 
AKARI. Thus, the shell detected with PACS is either a density 
enhancement due to a two-wind interaction or it could repre- 
sent the termin ation shock, si milar to the scenario proposed for 

Y CVn by [Libert et al.| ( |2057) . 

Another interesting finding with respect to Class III objects 
are the objects that have both a detached shell and a bow shock 
interaction region, with R Scl being the showcase object. Seeing 
the bow shock far out off the detached shell more or less excludes 
one of the explanations of detached shel ls, namely their in terpre- 
tation as being the wind-ISM interface ( Libert et al.|2007 ). 





HI CNM&WNM [DL90] (Eq. 4) 
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Fig. 11. ISM density derived from the observed stand-off dis- 
tance Rq as a function of the height above the Galactic plane, 
z. Different density scaling height models are also shown (e.g. 
Eqs . g] a nd|5). References : [DL9 0] : |Dickey & Lockman| ([1990] ) ; 
[DS94]: [Diplas & Savage| ( [T994l ); [L93]: |Loup et al.| ( |1993| ). 



7.4. Detached spherical shells 

The observed detached shell objects (Class III) show quite 
a range of expansion ages calculated out of their measured 



7.5. Comparison between predicted and observed Ro: 
Implications for the ISM density 

For a number of objects in our sample, previous studies have 
reported on the detection of wind-ISM bow shocks with IRAS, 
Spitzer or AKARI. The measured deprojected angular stand-off 
distance obtained with PACS agree very well with earlier re- 
sults (Table [7]). The predicted values for R (as function of ^Jn^) 
give ISM densities that are somewhat high er than those com- 
puted from Eq.[5] For a Ori, |Ueta et al. (2008 ) predict nu = 1.5 - 
1 .9 cm -3 which is indeed in line with nu = 1.9 cm" 3 derived from 
Eq. [5] but higher than th at inferred from th e stand-off distance, 
n n = 4.6 cm" 3 (Table |7}. |Ladjal et al.| ( |2oTo| ) derive n H >2 cm" 3 
for CW Leo, which is a factor of 2 higher than derived from both 
our predicted stand-off distance as well as from Eq.|5](^H = 0-5 
-1.0 cm" 3 ). 
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Fig. 12. Predicted stand-off distance as a function of distance 
(top) and as a function of height above the Galactic plane, z (bot- 
tom). All in units of parsec. Classes: I (stars), II (squares), III 
(circles), IV (diamonds), and X (crosses & triangles). 



the observed de-projected Rq with the theoretical Rq, taking nu 
as the unknown. The comparison of nu derived from Eq.[5]with 
that derived from the observed and predicted Rq shows that these 
generally agree within a factor of three (Fig. 11 >. In a few cases, 
the derived ISM density is an order of magnitude higher (n H « 5 
to 35 cm -3 ) than given by Eqs.|4]and[5] Possibly, for these cases, 
the star is moving through a diffuse cold medium (CNM) with 
typical densities ur ~ 10-100 cm -3 , an order of magnitude 
higher than those of the WNM (Sect. [5]). Alternatively, this ap- 
parent discrepancy could also be resolved if the ISM itself has a 
peculiar flow velocity of v*-vism with respect to the local stan- 
dard of rest. If the peculiar velocity between the two media is 
higher, the density required to arrive at the same stand-off dis- 
tance is lower, since v* oc 1 / y/njj. On the other hand, the de- 
rived ISM density also scales linearly with the adopted M, thus 
an order of magnitude over-estimate of the mass-loss rate leads 
to the same order of magnitude over-estimates of nu- This may 
explain the high nu for some objects with high adopted M, such 
as W Aql, but not easily for others, such as R Leo, which have 
already low values for M. For RLeo, the turbulent features in 
the bow shock indicate a higher space velocity which is at odds 
with the observed v* = 15 kms -1 . If indeed the relative star- 
ISM space velocity for R Leo is higher, the inferred ISM density 
will be correspondingly lower (now nu oc v~ 2 ). For example, 
for a much higher peculiar velocity of 45 kms -1 (more in line 
with the bow shock morphology), nu = 1.6 cm -3 which is close 
to % = 1.1 cm -3 given by Eq. [5] In Sect. [5] we roughly es- 
timated the LSR velocity of the local ISM and found that the 
corrections to v* are small (< 20%) for most objects, with a few 
exceptions. Only for two cases do these higher relative veloci- 
ties lead to much smaller values for the ISM density; nu = 10.1 
and 13.5 cm -3 for //Cep and UCam, respectively. The ISM ve- 
locity correction does not alter the high densities obtained for 
e.g. RLeo, WAql and EPAqr. Another explaination could be 
that this 'denser' medium is a remnant of an earlier, slower stel- 
lar wind. How and if the (tentative) binarity of both W Aql and 
EP Aqr plays a role is as of yet not clear. 



Adopting the values for the different parameters given in 
Tables [T] to [3] yields Ro less than 8' for all objects with detected 
bow shock interaction. For all detected wind-ISM interaction ob- 
jects, the predicted stand-off distances range from about 0.02 to 
0.9 pc, in line with the observed de-projected values. This gives 
some credibility in the predicted stand-off distances (from Eq.[2]) 
for the objects in Class X ( "non-detections "). These results sug- 
gest that additional bow shocks should have been detectable for 
other stars within 300 pc, in particular for R Dor, Y Lyn, RS Cnc, 
and possibly HD 100764 and RXBoo. We note that higher sen- 
sitivity observations could possibly still reveal the expected ex- 
tended emission. 

The uncertainty of the predicted stand-off distance depends 
on the different stellar and ISM parameters. The space velocity, 
v* (for nearby stars), as well as the wind velocity, v w are well de- 
termined observationally with small (-10%) uncertainties. The 
uncertainties in the mass-loss rates can be an order of magni- 
tude, in particular for the targets with inaccurate distance esti- 
mates (Sect. [5]) which introduces errors of a factor of three in R . 
Finally, the local ISM density is difficult to derive observation- 
ally. Without knowledge of the exact phase of the ISM the star is 
traversing, estimates of the density could be off by several orders 
of magnitude (see the discussion on densities, scale heights, and 
filling factors in Sect. [5]). Table [T] lists nu derived from equating 



7.6. Presence versus absence of wind-ISM bow shocks 

In the previous sections, we have established different morpho- 
logical classes of bow shocks and examined the basic physics 
giving rise to a bow shock and discussed the various parameters 
affecting their size and shape. In this section we explore differ- 
ent properties of the stellar objects as well as the ISM in order 
to understand which conditions are required to be able to detect 
bow shocks and/or detached shells around AGB stars and red 
supergiants. 

For 50 out of 78 (63%) AGB stars and red supergiant, we de- 
tect detached extended far-infrared emission suggestive of a bow 
shock, detached ring or irregular extended emission. Limiting 
the sample to nearby objects (d < 500 pc) eliminates many un- 
certainties and the "detection-rate" improves to 43 out of 56 ob- 
jects, or 78%. Restricting the distance further to, for example, 
d < 300 pc, only marginally improves the detection rate to 28 
out of 34 (80%). 

To test whether the "detection" and "non-detection" samples 
have similar or statistically different distributions as a function of 
d, z or nu, v*, i, M, and Rq, we use the two-sample Kolmogorov- 
Smirnov (K-S ) test. The results are given in Table [8] In order 
to eliminate the effect of large uncertainties in the distances, we 
perform the K-S test on a distance limited (d < 500 pc) sample. 
The K-S probabilities of "detections" versus "non-detections" 
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Table 8. Two-sample Kolmogorov-Smirnov test applied to dif- 
ferent stellar, circumstellar and interstellar parameters possibly 
relevant to "detection" and "non-detection" samples (distance 
d < 500 pc). 

Kolmogorov probability p a 





(IJIJV) 


(IJIJV) 


(III) 


(I-IV) 




vs. (Ill) 


vs. (X) 


vs. (X) 


vs. (X) 




n=(34,8) 


n=(34,16) 


n=(8,16) 


n=(42,16) 


critical value (a = 0.01) 


0.70 


0.56 


0.76 




d (pc) 


0.25 


0.18 


0.98 


0.28 


z(pc)/n H (cm- 3 ) 


0.99 


0.03 


0.19 


0.02 


v* (kms -1 ) 


0.43 


0.02 


0.09 


0.02 


inclination i (°) 


0.01 


0.52 


0.19 


0.68 


M(M Q yr 1 ) 


0.55 


0.48 


0.59 


0.57 


v w (kms -1 ) 


78 


0.74 


0.84 


0.78 


R predicted (pc) 


0.40 


0.07 


0.84 


0.09 


R predicted (arcmin) 


0.94 


0.02 


0.19 


0.02 



Notes. (a) A value close to unity indicates that the two samples (detec- 
tion and non-detection) have a high probability to be from the same 
parent distribution. A low probability suggests their distributions differ 
significantly, thus indicating the parameter is decisive in observing bow 
shock interaction in our survey The second to fifth column give the K-S 
results between different sub-sets of classes (see Table 4 for the classi- 
fication). The number of data-points in each (sub) set are given in the 
fourth row. 



are particularly low (and thus indicative that the two samples 
have different distributions) for nu/z, v*, and Ro. The corre- 
sponding distribution histograms for the most distinctive pa- 
rameters, riH, v*, and Ro, are shown for the different classes 

to 



(detection/non-detection) in Figs. [13] to [T5] This points towards 
the scenario that the presence of bow shocks is strongly depen- 
dent on the stellar velocity (relative to the local medium), the 
local ISM density and the resulting standoff distance. The latter 
sets the size of the bow shock region and is apparently deter- 
mined predominantly by v* and n^ and not so much by the star's 
mass-loss properties (cq. evolutionary phase). 

All objects with observed wind-ISM interaction zones 
(Class I to IV) have Ro < 1 pc, while many class X sources 
have Ro > 1 pc. The distribution of the observed de-projected 
Ro values is very similar to the distribution of predicted Ro for 
sources in Class I to IV. However, a quantitative comparison be- 
tween the predicted and observed stand-off distances does not 
reveal a strong correlation. 

Irrespective of the absolute distance (and thus apparent an- 
gular stand-off distance) it appears that a necessary, but not suf- 
ficient, condition for the detection of bow shocks and rings is 
that their physical size has to be smaller than about 1 pc. This 
points towards a physical effect (e.g. reduced surface brightness) 
instead of an observational limit, although it does not explain 
the absence of bow shock emission around nine of the observed 
nearby AGB stars (with known v* and Ro < lpc) included in 
Class X. 

Concluding, we find that bow shocks are detected and pre- 
dicted to occur for most nearby objects (< 500 pc) as well as 
for objects whose stellar and local ISM properties yield rela- 
tively small stand-off distances (< 1 pc). For more distant ob- 
jects the detection is (likely) hampered by lower sensitivity and 
lower spatial resolution, though we could not find any reasons to 
suggest that these objects would not have bow shock interaction. 
Indeed, several distant objects such as AQ And (-800 pc), S Cas 
(-940 pc), RZ Sgr (-730 pc) reveal evidence for wind-ISM in- 
teraction, but these detections might have been fortuitous. 




Fig. 13. Histogram of the "predicted" ISM density, nu, for the 
"detected" Class I, II, IV (top panel) and Class III (middle panel) 
"non-detection" Class X (bottom panel) objects. 




Class I, II, l\T 



Fig. 14. Histogram of the space velocity, v*, for the "detected' 
Class I, II, IV (top panel) and Class III (middle panel) "non- 
detection" Class X (bottom panel) objects. 




Fig. 15. Histogram of the standoff distance, Ro (pc), for the "de- 
tected" Class I, II, IV (top panel) and Class III (middle panel) 
"non-detection" Class X (bottom panel) objects. 



8. Conclusions 

This paper presents a morphological inventory of "wind-ISM" 
bow shocks and "wind- wind" interactions detected in the far- 
infrared with Herschel/PACS. Five main classes are identified: 

I "Fermata ". These objects are characterised by a large arc or 
shell-like structure spanning an angle of at least 90° . Several 
objects in this class show the presence of turbulence, i.e. 
Rayleigh-Taylor or Kelvin-Helmholtz instabilities. Some ob- 
jects, such as XPav, EP Aqr, and XHer, reveal the presence 
of back flowing material in the wake (behind) the star, lead- 
ing to a distinct bullet-shape. 
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II "Eyes". These are characterised by two arcs on opposing 
sides of the central object. These arcs are elliptical and 
non-concentric. In two cases the arcs connect and form a 
"lemon " shape morphology. 

III "Rings". This class includes the well known detached shell 
objects, typical of wind- wind interaction and dust pile-up. 
Younger, smaller rings have a co-spatial counterpart molec- 
ular gas ring, while the larger dust shells do not. 

IV "Irregular" . This class includes all sources which show ex- 
tended (non-detached), irregular emission at 70 or 160 fim. 

X "Non-detection" . All objects for which no large extended 
far-infrared emission has been observed (c.q. resolved) in the 
form of a bow shock or ring with the PACS instrument. 

Oxygen rich stars give rise predominantly to "fermata" or 
"irregular" morphologies. There is tentative evidence that all 
"ring" stars are thermally -pulsing AGB stars. We identify a few 
cases for which both a detached (inner) shell is found within the 
(outer) bow shock, suggesting that the detached shell is not due 
to the wind-ISM interaction, but suggest rather a wind- wind sce- 
nario. 

The presence or absence of bow shocks is determined by the 
stand-off distance which depends on relative values of the stel- 
lar parameters (v*, v w and M) as well as those of the ISM (^ H ) 
as given by the theoretical approximation as well as the hydro- 
dynamical simulations. The distribution of stand-off distances 
for the entire survey sample shows a clear separation between 
objects with and without detected extended emission. Indeed 
most of the objects (44) assigned to Classes I to IV have Ro < 
1.0 pc. No extended bow shock or detached shell emission has 
been detected for the 15 objects with predicted Ro > 1 pc and/or 
d > 1000 pc. Limiting the sample to d < 500 pc, we derive a 
Ro < 1.0 pc, and thus predict the presence of bow shocks, for 8 
stars in the "non-detection" class: TWHor, VEri, RDor, RLep, 
RS Cnc, HD 100764, RYDra, and RXBoo. Many of these show 
extended shells at about 1 MJy sr" 1 sensitivity in Spitzer and 
AKARI maps (private communication: T. Ueta for Spitzer and 
H. Izumiura for AKARI). 

The angular size of the predicted stand-off distances for 6 
other objects in Class X are larger than the obtained PACS im- 
age maps. In particular, all-sky survey missions, such as AKARI 
( [Murakami gjjggggT) and WISE ( [Wright et al.1[20T0| ) could 
help in observing extended bow shocks and shells around these 
objects. Observations of bow shocks and detached shells for 
distant (>1 kpc) AGB stars and red supergiants require thus 
higher sensitivities to detect extended emission, but not neces- 
sarily higher spatial resolution as most have predicted stand-off 
distances larger than V . 

The observed infrared emission indicates the presence of 
moderate amounts of dust (2 x 10" 7 to 2.5 x 10" 4 M ), implying 
only up to a few tenths of solar masses of dust and gas in the 
bow shock interaction region. The observed dust and gas masses 
are similar to the potentially swept-up ISM material. 

This survey represents only a first step in fully characteris- 
ing and understanding the formation and shaping of bow shocks 
and detached shells around AGB stars and red supergiants. It is 
clear that additional simulations trying to represent as closely 
as possible the different observed morphologies are required in 
particular to understand the formation of instabilities. Detailed 
observations of the bow shock spectral energy distribution (dust) 
and line emission (gas) is needed to quantify the physical condi- 
tions and composition of the material in these shocked regions. 
Potentially, detailed observations of bow shocks around larger 
samples of AGB stars can be used to independently probe the 



local ISM density (and possibly the magnetic field) distribution 
in the Galaxy. 
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